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Abstract 

We consider the problem of enumerating permutations in the symmetric group on n elements 
which avoid a given set of consecutive pattern S, and in particular computing asymptotics as n 
tends to infinity. We develop a general method which solves this enumeration problem using the 
spectral theory of integral operators on L 2 ([0, l]™ 1 ), where the patterns in S has length m + 1. 
Krem and Rutman's generalization of the Perron-Frobenius theory of non-negative matrices plays 
a central role. Our methods give detailed asymptotic expansions and allow for explicit computation 
of leading terms in many cases. As a corollary to our results, we settle a conjecture of Warlimont 
on asymptotics for the number of permutations avoiding a consecutive pattern. 

1 Introduction 

In this paper, we study integral operators of the form 

(Tf)(xi,...,x m )= x(t,xx, . . . ,x m ) • f(t,xi, . . . ,x m -i)dt (1.1) 
Jo 

and their application to enumerating permutations that avoid a consecutive pattern. Here x is a real- 
valued function on [0, l] m which takes values in [0, 1] and is continuous away from a set of measure 
zero in [0, l] m+1 . As we will show, operators of this type arise naturally when counting permutations 
that avoid a consecutive pattern of length m + 1. 

To define the enumeration problem, let & n denote the symmetric group on tl elements. For 7r G &n 
we write ir = (ttitt2 ■ ■ ■ ir n ) where the tt^ are the integers from 1 to n. For x = (xi, . . . , x n ) G W 1 with 
Xi ^ Xj for all distinct i and j, we denote by n(x) the unique permutation ir G & n with tt-i < 7Tj if and 
only if Xi < Xj. A pattern of length k is an element a of If it is a permutation of length n > k, it 
avoids a if n(7r J ', itj+i, ■ ■ ■ , Kj+k-i) 7^ & for all j with \ < j < n — k + 1. More generally, if S C <& k we 
say that 7r avoids S if tt avoids each a G S. That is, S is the set of forbidden patterns. 

For a subset S of <5 m +i, denote by a n (S) the number of permutations ir £ & n that avoid S. 
Observe that for n < m we have a n (S) = n\. Our goal is to compute asymptotics of a n (S) as n tends 
to infinity. Throughout the paper we will assume that m > 2. 

For S C 6 m+ i we define a function \ on [0, l] m+1 by setting x( x ) = if Xj = Xj for some distinct 
indices i and j, and otherwise setting 

= { i n(x) t s, 
x[x> \ o n(x) g s. { > 

Finally, let T be the operator of the form (II. ip . The standard inner product on L 2 ([0, l] m ), is defined 

by 



(f,g) = / f(x) -g(x)dx, 

'[0,l] m 
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and hence the I? norm is given by ||/|| = -J (f, /). Moreover, let 1 denote the constant function 1 
on [0, l] m . It is straightforward to prove (see Proposition 12, 10] ) that the formula 

= (T n - m (l),l) (1.3) 

n! 

holds for n > m. Note that the left-hand side is the probability of selecting a permutation ir S & n at 
random that avoids the set S. 

The asymptotic behavior of powers of a bounded linear operator is determined by its spectrum. 
Recall that if A is a bounded linear operator from a Hilbert space to itself, the resolvent set of A is the 
set of all z S C so that (zl — A) -1 is also a bounded operator. The complement in C of the resolvent 
set is the spectrum of A, denoted cr(A), and the spectral radius of A is given by 

r(^) = sup{|A| :\£a(A)}. 

The peripheral spectrum of A is the intersection of cr(A) and the circle of radius r{A) in the complex 
plane. As we will see, the peripheral spectrum of the operator T consists at least of a real eigenvalue 
at r(A), although this need not be the only eigenvalue in the peripheral spectrum. Finally, define the 
adjoint of an operator A to be the operator A* that satisfies (/, A*(g)) = (A(f),g). 
Using spectral theory, we obtain: 

Theorem 1.1. Let S be a set of forbidden patterns in & m+ i. Then the nonzero spectrum of the 
associated operator T consists of discrete eigenvalues of finite multiplicity which may accumulate only 
at 0. Furthermore, let r be a positive real number such that there is no eigenvalue ofT with modulus r 
and let Ai, . . . , A& be the eigenvalues ofT greater in modulus than r. Assume that Ai, . . . , A& are simple 
eigenvalues, with associated eigenf unctions (fi and that the adjoint operator T* has eigenf unctions ipi 
corresponding the eigenvalues Aj. Then we have the expansion 

a n {S)/n\ = (T"~ m (l), 1) = V to" 1 )' Mi) . X n m + 0{r ny ^ 

Moreover, when the operator T has a positive spectral radius, that is, r{T) > then the spectral radius 
is an eigenvalue of the operator T. 

Observe that when T has spectral radius 0, this result is not useful. 

In applications the error term in equation f|l .4[) can be selected to be the modulus of the next 
eigenvalue, that is, 

k 

a n {S)/n\ =J2 c i- A i + °(\ X k+i\ n ), (1-5) 
i=i 

where the eigenvalues satisfy |Ai| > • • • > |A&| > |Ajt+i| and A^+i is a simple eigenvalue. This follows 
from letting r be smaller than |Afc+i| given a few more terms in equation (II. 4j) and increasing the error 
term. 

For certain sets of patterns, we can show that there is a unique largest eigenvalue with respect 
to modulus and that this eigenvalue is simple and positive. Thus we obtain the leading term of the 
asymptotics of a n (S)/n\. A sufficient condition for the peripheral spectrum of T to consist only of 
the simple eigenvalue r(T) is as follows. If (X, fi) is a measure space and / E L 2 (X, /i) is a real-valued 
function, we will say that / > if f{x) > for almost every x £ X, and / > if f(x) > for almost 
every x. A bounded operator A on L 2 (X,/j.) is positivity improving if for any / > different from 
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there is an integer k (possibly depending on /) so that A k f > 0. Krein and Rutman |15[ Theorem 6.3] 
showed that, for such an operator, the spectral radius r(A) is a simple eigenvalue r(A) and all other 
eigenvalues A are smaller than r(A), that is, |A| < r{A). Furthermore, the associated eigenfunction ip 
is positive almost everywhere. Finally, the adjoint operator A* also has a largest real, positive and 
simple eigenvalue at r(A) and an almost everywhere positive eigenfunction ip. We summarize this 
discussion in: 

Theorem 1.2. If the operator A is positivity improving then its largest eigenvalue is real, positive and 
simple. 

A sufficient condition for T to be positivity improving can be formulated in combinatorial terms as 
follows. If 5 is a set of patterns, let Gs be the graph with vertex set & m and a directed edge from tt to 
a if there is a permutation r £ 6 m+ i\5 with n(ri, . . . , r m ) = tt and II(t2, . . . , r m+ i) = a. The graph 
G$ is known as the graph of overlapping permutations. Recall that a graph G is strongly connected if 
any vertex of G is connected to any other vertex by a directed path. 

Theorem 1.3. Suppose that the graph Gs is strongly connected and that the monotone permutations 
12 • • • (m + 1) and (m + 1) • • • 21 do not belong to S. Then the operator T is positivity improving. 

As a corollary we have the following result, proving a conjecture of Warlimont [22]. See also 
Theorem 4.1 in [9]. 

Corollary 1.4. Let S consists of a single permutation a. Then the asymptotics of a n (S) is given by 
a n (S)/n\ = c • A™ + 0(r n ), where c, A and r are positive constants such that A > r. 

Proof. When the permutation a differs from the two monotone permutations 12 • • • {m + 1) and (m + 
1) • • • 21, Theorem 11.31 applies to the forbidden set S = {cr}. The two remaining cases are equivalent, 
and can be settled by Theorem 11.71 □ 

Another application of Theorem 11.31 is the following. Call a permutation tt in © n decomposable 
if there exists an index i such that 1 < i < n — 1 and 7r(l),...,7r(i) < i (equivalently i + 1 < 7r(i + 
1), . . . ,7r(n)). A permutation that is not decomposable is called indecomposable. These permutations 
are also known under the terms connected and irreducible. 

Theorem 1.5. Let S be a subset of G m+ i such that each permutation in S is indecomposable and the 
monotone permutations (m + 1) • • • 21 do not belong to S. Then a n (S) has the asymptotic expression 
a n (S)/n\ = c • A™ + 0(r n ), where c, A and r are positive constants such that A > r. 

More generally we can characterize the spectrum of T in terms of another graph associated to S. 
To define it, let A n denote the set of points x = (xi,...,x m ) E (0, l) m with Xi ^ xj for i ^ j 
and Il(a;) = tt. The graph Hs has vertex set U 7T ^Q m A n and directed edges from (xi,...,x m ) to 
(x2, • • • ,x m+ i) if x\ ^ x m+ i and n(xi, . . . ,x m+ i) ^ S. We show that if the graph Hs is strongly 
connected there is an upper bound on the length of the directed path connecting any two vertices. 
We define the period of a strongly connected graph G as follows. Fix a vertex v of G and, for k a 
non-negative integer, let be the set of all vertices in G that can be reached from v in exactly k 
steps. The set Q of all k with v G is a semigroup and generates a subgroup dZ of Z. The integer 
d is the period of the graph G. Note that, if G is strongly connected, then G has period d for some 
positive integer d. Finally, a graph is called ergodic if it is strongly connected and has period 1. 
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Theorem 1.6. Suppose that S is a set of forbidden patterns and that the graph Hg is strongly con- 
nected with period d. Then the operator T has positive spectral radius r(T) and T has a simple 
eigenvalue A = r(T) with strictly positive eigenf unction. Moreover, the spectrum of T is invariant 
under multiplication by exp(2 • tt ■ i/d). 

In case when the period is 1, we have the conclusion: 

Theorem 1.7. Suppose that S is a set of forbidden patterns such that the graph Hs is ergodic. Then 
the operator T is positivity improving. That is, the operator has a unique largest eigenvalue which is 
simple, real and positive and the associated eigenfunction is positive. 

For certain explicit patterns, we can compute the spectrum and eigenfunctions of T and obtain 
sharp asymptotic formulas for a n (S)/n\. 

Example 1.8. When S is empty, directly we have ct n (S) = n\ for all n > 0. The associated operator 
only has one non-zero eigenvalue, namely 1, with eigenfunction and adjoint eigenfunction 1. This is 
the only case we know where the number of eigenvalues is finite. 

Example 1.9. If S = {123} we show that the operator T has a trivial kernel and spectrum given by 
i x k} k& where 

A ^ 
k 2.*.(* + i)- 

Furthermore all the eigenvalues are simple. We also compute the eigenfunctions of T and the adjoint 
operator T* and obtain 

a„(123) ( 1 \ yll+1 



nl «P^J-*r+°(|A-lH CM) 

where Ao = r(T) and A_i is the next largest eigenvalue in modulus. For more terms in the asymptotic 
expansion see Theorem 15.41 

Example 1.10. If S = {213}, we show that the nonzero eigenvalues of the operator T are the roots 
of the equation 

erf ivh) = V? 

which has the unique real root Ao = 0.7839769312 .... Moreover, Ao is the largest root in modulus of 
the equation. We then have 



a„(213) 



exp 



n! ~' r V 2 • Xl 



Xq +1 + O (|Ai| n ) (1.7) 



where Ai^ = 0.2141426360 . . . ± 0.2085807022 . . . • i are the next two largest roots of the eigenvalue 
equation. See Section [6] for the calculations. 

Example 1.11. If S = {123, 321}, the numbers a n (S) are given by a n (S) = 2E n for n > 2 where E n 
is the nth Euler number. In this case, we can use spectral methods to obtain the classical convergent 
expansion 

§=2. y, (^y 1 - <L8 » 

j odd 
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This formula was derived by Ehrenborg, Levin, and Readdy [8j (Corollary 4.2) by using Fourier series. 
In this case the spectrum of the operator T real and invariant under the reflection A i— > —A; in particular 
the two largest eigenvalues are ±2/tt. Moreover, the matrix U introduced in Proposition ^. Ill is similar 

to the matrix ^ for a cyclic permutation of order two. 

Example 1.12. Let S = {123, 231, 312}. In this case the asymptotic expansion converges and we can 
conclude that the number of such permutations in 6 n is given by n ■ E n _i, a result due to Kitaev and 
Mansour [13]. See Section [71 

It is easy to find examples of patterns S for which p{T) = 0. 

Example 1.13. Let S = {132,231}. An S'-avoiding permutation has no peaks (viewed as the graph 
of a function from {1, . . . , n} to itself) and it is easy to see that a n (S) = 2 n ~ l . However, observe that 
the proportion a n (S)/n\ is subexponential. It is straightforward to verify that the operator T has no 
non-zero eigenvalues. Also note that the graph Gs is not strongly connected. 

Example 1.14. Let S = {123,213,231,321}. The directed graph Gs is strongly connected but the 
monotone permutations 123 and 321 are excluded. In this case a n (S) = 2 for all n > 2. 

We close our introduction by a brief overview on the subject of pattern avoidance in permutations 
(for more details we refer to [3J). The "classical" definition of a pattern is slightly different than one 
provided above. We say that a permutation n avoids a pattern a if tt does not contain a subsequence 
which is order-isomorphic to a. The study of such patterns originated in theoretical computer science 
by Knuth [14]. However, the first systematic study was done by Simon and Schmidt [18] . who com- 
pletely classified the avoidance of patterns of length three. Since then several hundred papers related 
to the field have been published. 

One of the most important results in the subject is the proof by Marcus and Tardos [16] of the 
so-called Stanley-Wilf conjecture related to the asymptotic behavior of the number of permutations 
that avoid a given pattern. It states that for any permutation a there exists a constant c (depending 
on a) such that the number of the permutations of length n that avoid a is less than c n . 

In this paper we also study asymptotic behavior of permutations avoiding patterns, but we con- 
sider consecutive patterns, occurrences of which correspond to (contiguous) factors, rather than sub- 
sequences, anywhere in permutations. Simultaneous avoidance of consecutive patterns of length 3 
is studied in [12] using direct combinatorial arguments. Other approaches to study consecutive pat- 
terns were introduced recently which include considering increasing binary trees [10] by Elizalde and 
Noy, symmetric functions [T7] by Mendes and Remmel, homological algebra [5] by Dotsenko and 
Knoroshkin, and graphs of pattern overlaps [T] by Avgustinovich and Kitaev. 

In [10] asymptotics for the following consecutive patterns is given: 123, 132, 1342, 1234 and 1243. 
These results are obtained by representation of permutations as increasing binary trees, then using 
symbolic methods followed by solving certain linear differential equations with polynomial coefficients 
to get corresponding exponential generating functions, and, finally, using the following result (see [11[ 
Theorem IV. 7] for a discussion). 

Theorem 1.15. If f(z) is analytic at and R is the modulus of a singularity nearest to the origin in 
the sense that 

R = sup{r > : / is analytic in \z\ < r} , 
then the coefficient f n = [z n ]f(z) satisfies 

limsup|/ n | 1 / n = l/i?. 

n— >oo 
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This differ from our method, as we obtain detailed asymptotic expansions that allows for explicit 
computation of leading coefficients in many cases. As special cases of our results, we get more detailed 
asymptotics for some of the results of Elizalde and Noy [10] . Also note that the associated generating 
function is related to our operator T by the identity 



MS) ■ Z -, = 1 + • • • + z m ~ l + z m • ((/ - zT)-\l), 1) . 

n>0 



From this identity it is clear that the radius of convergence of the generating function is determined 
by the spectrum of T. 

The underlying principle that makes our method work is that one can pick a permutation in & n at 
random with uniform distribution, by picking a point (x\, . . . , x n ) in the unit cube [0, l] n and applying 
the function II. This method has been used in [8] to obtain quadratic inequalities for the descent set 
statistics and in [7] to enumerate alternating 2 by n arrays. 

2 The operator T 

We now begin our study of the operator T. 

2.1 Connection with pattern avoidance 

Recall that S is a collection of forbidden patterns of length m+ 1, that is, S is a subset of <5 m +i. The 
function x is defined on the unit cube [0, l] m+1 by 

/ 1 U(x) $ S, 
nx) \ n(a;) G s, 

and the operator T on L 2 ([0, l] m ) by 

(Tf)(xi,...,x m )= / x(t,xx,...,x m ) ■ f(t,xi,...,x m -x)dt. 
Jo 

For n > m, define Xn on the n-dimensional cube [0, l] n by 

n—m 

Xn(xi, • • • ,x n ) = Y[x(xj,---,x m+j ). (2.1) 

3=1 

This allows us to express powers of our operator T. If k < m then 

(T k f) (x) = / Xm+k{h, ■ ■ ■ ,t k ,xi, . . . ,x m ) ■ /(ii, ...,t k ,xt,... ,x m - k )dti ■■■dt k 
v ' ifo.il* 



(2.2) 



while if k > m, 



[T k f)(x)= / Xm+k(h, ■ ■ ■ ,tk,xi, . . . ,x m ) ■ f(ti, . . . ,t m )dt\ ■ ■ ■ dt k (2.3) 

v ' i[o,i] fc 

Proposition 2.1. The number of S-avoiding permutations in G n for n > m is given by a n (S) = 
n\ ■ (T n " m (l),l). 
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Proof. Equation (I2.ip allows us to conclude that 

Xnix) = 

Hence by integrating over the n-dimensional cube 



1 if II(x) avoids S, 
if II(x) does not avoid S. 



[0,1]" 



a n (S) 

Xn(x)dx - 



ni 



since each simplex in the standard triangulation of [0, 1]" corresponds to a permutation tt S <5 n and 
each such simplex has volume 1/n!. Using this observation and the identities (12. 2p and (12. 3p . we can 
rewrite the integral as (T n_m (l), 1). □ 

Lemma 2.2. The operator T is a bounded operator with norm at most 1. 

Proof. Apply the Cauchy-Schwarz inequality to equation in the variable t to obtain 

1 2 

X(t, xi,..., x m ) ■ f(t, xi,..., x m -i)dt 



o 



< / \x{t,xi,...,x m )\ dt- / |/(t,xi, . . . , x m _i)\ dt 
Jo Jo 

< / \f(t,xi, . . . ,x m -i)\ 2 dt. 



Now integrating over the variables x\, . . . , x m , we obtain that ||T(/)|| 2 < ||/|| 2 proving the bound. □ 
Lemma 2.3. The adjoint operator T* is given by 

(T*f)(x)= / x(xi,---,x m ,u)f(x2,...,x m ,u)du. 
Jo 



Proof. Since x is real- valued, we have that 

{T(f),g) = / x(t,%i, ■ ■ ■ ,Xm) ■ f(t,xi, ■ ■ ■ ,x m -x)dt ■ g(xi,... ,x m )dxx ■ ■ -dx r 

J[o,i] m Jo 



fit 

, X\ , . . . , X m —\ )dt- I x(t,xi,...,x m )-g(x 1 ,...,x m )dx m dtdx 1 ---dx m ^i, 
'[o,i] m Jo 

proving the lemma. □ 

The spectrum of an adjoint operator A* is given by conjugate of the spectrum of A. That is, if 
A is an eigenvalue of A, then A is an eigenvalue of A*. However, for our operator T, the complex 
eigenvalues come in conjugate pairs. This is proved using that x is real- valued. 



Lemma 2.4. We have that T(f) = T(f) and T*{f) = T*(f). Hence if X is an eigenvalue of T with 
eigenfunction ip, then A is also eigenvalue of T with eigenfunction Tp. Similarly, if A is an eigenvalue 
of T* with eigenfunction ip, then A is also eigenvalue of T* with eigenfunction ip. 
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2.2 Eigenvalues and asymptotic expansion 



Since the operator T is bounded, we can use spectral analysis in order to explore the operator T. We 
refer the reader to Dunford and Schwarz [6, Chapter VII] for a more detailed exposition. 

As we defined in the introduction, the resolvent set s(T) is the set of complex numbers z such the 
operator (zl — T) _1 exists as a bounded operator. The spectrum o~(T) of T is the complement of the 
resolvent set ?(T). 

The index of a complex number A is the smallest non-negative integer v such that the equation 
(XI — T) u+1 f = implies (XI — T) u f = for all functions /. Informally speaking, for operators on 
a finite-dimensional vector spaces, the index of an eigenvalue is the size of the largest Jordan block 
associated with that eigenvalue. A point A in the spectrum is called a pole T of order v if the function 
R(z; T) has a pole at A of order v. Theorem 18 in [6, Section VII. 3] states that the order of a pole A 
is equal to its index. 

Define the operator E(X) by the integral 



where C is a positive oriented closed curve in the complex plane only containing the eigenvalue A from 
the spectrum cr(T). It follows from [61 Section VII. 3] that E(\) is a projection. 

Lemma 2.5. The operator T m is compact. 

Proof. The operator T m has the form 



Since X2m is a bounded function we conclude that T m is a Hilbert-Schmidt operator, and hence a 



Using Theorems 5 and 6 in [6l Section VII. 4] we conclude: 

Theorem 2.6. The spectrum of T is at most denumerable and has no point of accumulation in the 
complex plane except possibly 0. Every non-zero number X in a(T) is a pole ofT and has finite positive 
index. For such a number X the projection E(X) has a non-zero finite dimensional range and it is given 
by {/ 6 L 2 ([0, l] m ) : (XI - T) v f = 0} ; where v is the order of the pole X. 

Recall that an eigenvalue A is simple if the range of E(X) is one-dimensional. That is, the eigenvalue 
equation Xtp = Tip has a unique solution up to a scalar multiple and the generalized eigenvalue equation 
Xf = Tf + ip has no solution. 

Lemma 2.7. Let X be a simple eigenvalue of the operator T with associated eigenfunction (p. Let ip 
be the eigenfunction of the adjoint operator T with eigenvalue X. Then the projection E(X) is given by 



Proof. Since the eigenvalue A is simple, the range of the projection E(X) is one-dimensional and 
spanned by the eigenfunction (p. Since the projection is continuous we may assume that it has the 
form E(X)(f) = (f,a) ■ (p for some function a. It is straightforward to observe that since E(X) is a 
projection, that is, -E(A) 2 = E(X), we have that (<p,ac) = 1. 





compact operator. 



□ 
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Now the adjoint operator E(X)* is given by E(X)*(g) = (g, </?) • a since 

(g, E(X)(f)) = (g, (/, a)-<p) = (a, /) • (g, <p) = ((g, if) ■ a, /) . 

Hence a belongs to the range of the adjoint operator, that is, it is multiple of the eigenfunction ip. In 
fact, we observe that a is given by Up, ip) • tp. □ 

Proof of Theorem By analytic functional calculus we can evaluate the operator T n ~ m by inte- 
grating in the complex plane; see Theorem 6(c) in [61 Section VII. 3]. We have 



1 



2vri T\ z \=r zl — T 



dz, 



where R is greater than the spectral radius of T and we orient the circle in positive orientation. 

Let a be the set {Ai, . . . , A&} and let E(a) denotes the sum of the projections E{X\) + • • • + E(Xk)- 
By Theorem 22 in [6j Section VII. 3] and that the eigenvalues Ai, . . . , are simple, we have that 



T n - m ■ E(a) = J2 E(Xi) ■ A™"™ 



i=l 



We can estimate the operator T n rn ■ (I — E(a)) by shrinking the path of integration to a circle of 
radius r 

T ". (I-E{a)) = — ■ 4) — —(I:. 



2m 



We bound this integral by 

\\T n ~ m ■ {I - E{o))\\ 



1 




2tH ' 


J\z 



\z\=r 



< 

~ 2vr 



\z\=r 



z n- 


-m 


-_ r zl - 


- T 


z n- 


-m 


-. r zl - 




- T 


1 




zl - 


T 



dz 

dz ■ r n - m 



-i 



< sup (zl — T) 

\z\=r 

= 0{r n ), 

where the last equality follows from that the supremum does not depend on n. Hence the inner 
product (T n_m • (I — E(a))l, 1) is also bounded by 0(r n ). Thus we conclude that 

(T"- m l,l) = (T n - m E(a)l,l) + (T n ~ m - [I - E(a))l,l) 

k 



Y J {E{X l )l,l)-Xr m + 0[r n ) 

+ 0(r n ). 



i=l 
k 



i=l 



□ 
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If the eigenvalue A is not real then A is also an eigenvalue. The two terms corresponding to these 
two eigenvalues can be combined as follows. Let A = r • e 1 ' 8 and (ip, 1) • (l, ip) / up, tp) = s ■ e 1 '^ . Then 
we have 

A „_ m | (p,l)-(l,ip) yn- m = 2 .ft( a . e *'P . r n-m , e i-{n~m)-e 

(<p,W) (WiVO ^ 

= 2-s-r n - m -cos(f3 + (n-m)-6). 
2.3 Bounds on the norm and spectral radius 

Proposition 2.8. Let T be an operator of the form M.l\) and suppose that T has a positive spectral 
radius r(T). Then T has a positive eigenvalue A = r(T) with non-negative eigenfunction ip. Further- 
more, the eigenfunction ip is almost everywhere positive. Similarly, the adjoint operator T* also has 
A = r(T) as an eigenvalue with an almost everywhere positive eigenfunction ip. 

Proof. Let K be the cone of non-negative functions in L 2 ([0, l] m ). Since T preserves this cone and T 
has nonzero spectral radius, it follows from Theorem 6.1 of [15] and the fact that the cone of positive 
functions is self-dual that T and T* both have A = r(T) as an eigenvalue with at least one strictly 
positive eigenfunction. □ 

For an example of an eigenfunction ip corresponding to the largest eigenvalue, taking the value 
on a set of measure 0, see Proposition 17. 21 where 123, 231, 312-avoiding permutations are discussed. 



Lemma 2.9. For k > m the norm of the operator T k is bounded above by y/ a m +fc (S) / (m + k) ! . 

Proof. Applying the Cauchy-Schwarz inequality to equation (|2.3p in the variables t\ through tk, we 
obtain 

2 

Xm+k (t 1 , • • • , tk , X\ , . . . , X m ) • f(ti, . . . ,t m )dt\ 



[0,1] 



< / \Xm+k(tl,. ■ ■ ,t k ,xi, . . . ,x m )\ dt ■ \f(ti,...,t m )\dt 
J[0,l] k J[0A] k 



2 



Xm+k(h, ■ ■ ■ ,t k ,xi, . . . ,x m )dt ■ / \f(tx,...,t m )\ dt. 

[0,l] fc J[0,l] m 

Now integrating the variables x\,...,x m over [0, l] m we have ||T fc (/)|| 2 < a m+k {S)/(m + k)\ 
proving that ||T fc || < y/a m +k(S)/(m + A;)!. □ 

Proposition 2.10. For a nonempty set of forbidden patterns S we have that spectral radius of the 
operator T is less than 1, that is, r{T) < 1. 

Proof. We have a2 m (S) < (2 • m)\ implying that ||T m || < 1. The inequality ||T m ' n || < ||T m || n implies 
that r(T m ) = lim 

n — >oo \\T mn \\ ' n < 1 strictly. The result follows by taking the ?nth root. Q 

Proposition 2.11. Let T be an operator of the form M.l\) and suppose that p = r(T) > 0. Then the 
operator T admits a decomposition of the form 

T = p-U + W 

where UW = WU = 0, U has finite- dimensional range, U maps the interior of the cone of positive 
functions into itself, the operator W has spectral radius less than p, and the eigenvalues of U are roots 
of unity including 1 . 
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Proof. The orthogonal decomposition follows from Theorem 8.1 of [JS] applied to the operator A = 
p~ l -T. □ 

In particular, the leading behavior of powers T n is determined by the spectral radius r(T) and the 
finite-rank operator U. 

Theorem 2.12. Let S be a set of forbidden patterns. Then spectral radius of the operator T is given 
by 

r(T)= lim 

n — s>oo \ n! 



Proof. Suppose first that r(T) = 0. From the inequality |(T n l,l)| < ||T n || we immediately conclude 
that (a n (S')/n!) 1//n tends to as n goes to infinity. If r(T) > 0, then by Proposition 12.11] we have 

^ifl = r (T) n OJ n l, 1) + (W n l, 1) 
n! 

where the second term obeys the estimate 

\(W n l,l)\<(r(T)-e) n 

for some e > and all sufficiently large n. Moreover, (U n l, 1) is periodic in n and strictly positive 
since U maps the interior of the cone of positive functions into itself. Thus (U n l, 1) is both bounded 
above and below by strictly positive constants. It follows that lim, t — >OQ {U n l,l) 1/n = 1 so that 
lim n >ao (a n (5)/n!) 1//n = r(T) as claimed. □ 

This result extends Theorem 4.1 of [9], where consecutive patterns consisting of a single permuta- 
tion were considered. Moreover, Theorem 12.121 characterizes the limit in terms of a spectral quantity 
which can be computed in many cases of interest by solving the eigenvalue problem for the integral 
operator T. 



3 Associated graphs 
3.1 The directed graph Hs 

In this section we study the spectrum of T using the infinite graph Hs described in the introduction. 
Recall that A n denotes the open subset of (0, l) m with Xi ^ Xj for i 7^ j and Xi < Xj if and only if 
7r(«) < vr(j), and let 

Thus the complement of X consists of those points x with X{ = Xj for at least one pair of distinct 
indices i and j and hence is a set of measure zero. The graph Hs has vertex set X. Recall that the 
directed edges of Hs connected points x and y in X with Xj + \ = yj for 1 < j < m — 1, x\ ^ y m , 
and II(:ei, . . . , x m , y m ) £ S. It follows from the definition that x(t, x±, . . . , x m ) = 1 if and only if there 
is a directed edge from (t, x±, . . . , x m -i) to (xi, . . . , x m ). That is, the function \ encodes the edge 
information of the graph Hs ■ 

The next lemma connects the graph Hs to mapping properties of the operator T. It will be used 
to show that the operator is positivity improving. 
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Lemma 3.1. Suppose that x,y G X and that there is a directed path from x to y of length k > m. 
Suppose further that f is a non-negative continuous function such that f is non-zero in a neighborhood 
ofx. Then (T k f)(y) > 0. 

Proof. Assume that the directed path is 

x = (xi, . . -,x m ) — ► (x 2 , ■ ■ -,x m+ i) — > ■ ■ • — ► (x k+1 , . . .,x k+m ) = y. 

Let e be the minimum of the following finite set 

{\xi — Xj\ : 1 < i < j < k + m, j — i < m} U {xj, 1— Xi : 1 < i < k + m} . 

Observe that e > by the definition of X. Let 5 = e/3. For Sj G [xi — S, x^ + <5], 1 < i < k + m, we 
have that 

(si,...,s m ) — > (s 2) ... 

is also a directed path in Hs- It follows that Xk+m(si, ■ ■ ■ , Sk+ m ) = 1 for all such s. Using (|2.3p we 
may estimate 

"xi+6 rxk+S 



. . rxi+6 rx k +6 

(T k f) (y) > / ••• / f(t l ,...,t m )dt 1 ---dt k 



= (2S) k - m / ••• / f(t 1 ,...,t m )dt 1 ---dt m 

Jxi—5 Jxm — S 

> o, 

where in the last step we have used the positivity of / in & neighborhood, of (3^1, ■ ■ • ; %m 
Proposition 3.2. Suppose that Hs is strongly connected with period d. Then there is a decomposition 

d-l 



X=[jY i 



i=0 

of X into disjoint sets Y{ with the property that T : L 2 (Yi) — > L 2 (Yi + \), where Yd = Yq. 

Proof. Pick a base vertex v of Hs- Let X k be the set of all vertices in Hs that can be reached from v 
in k steps, and let Q be the subset of the non- negative integers defined by 

Q = {k : v£X k }. 

Then Q is a semigroup under addition and generates a subgroup of the integers Z. A subgroup of Z 
has the form dL for some positive integer d; in this case, d is the period of the graph Hs- 
Now define 

Y t = |J Xj. 

j: j=i mod d 

Observe that every directed edge in the graph Hs goes from some Yi to the next Y,- L+ i (with addition 
modulo d). Also, observe that the sets Y{ are pairwise disjoint. 

We claim that each Yi is open. To see this, suppose that y E Y^. Pick a path from some vertex x 
to the vertex y having length greater than 2m. We can perturb this path in a small neighborhood of 1 
y using a variant of the argument used at the beginning of Lemma IBTTl and conclude that Yi is open. 
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Now pick a permutation it € & m - Since the sets Yq,..., l^-i are pairwise disjoint we have that 
A n = A n n X = U^~ 1 (A 7r (~)Yi). Note A^ is a connected set and A n n 3^ are all open. A connected 
set can only be the disjoint union of one open set and hence there exists a unique index i such that 
Ajr C Y{. Hence each set Yi is the disjoint union of the sets A n . 

Finally, suppose that / is a continuous function on [0, l] m with support in the set Y^. We claim 
that Tf is supported in the next set Yi + ±. To see this, note that x(i,xi, . . . ,x m ) = 1 if and only if 
there is a directed edge from (t, x±, . . . , x m _i) to [x\, . . . , x m ). Assuming that f(t, x\, . . . , x m -i) is 
non-zero where (t, xi, • • • , x m —i ) belongs to the set Yi. Since / is continuous function, / is supported 
in a neighborhood of the point (t, xi, . . . , x m _i). By applying the definition of the operator T, we have 
the function Tf supported in a neighborhood of the point (xi, . . . , x m ) in the next set Yi + \. That is, 
Tf is supported in Yi + \. □ 

The next lemma is a straightforward consequence of the definition of the graph H$ and hence the 
proof is omitted. 

Lemma 3.3. Suppose that a : (0, 1) — > (0, 1) is a strictly increasing function. If 

(xi, . . . , X rn ) y (X2, • • • , Xm+l) > • • • y (Xk+l, • • • , Xk+rn) 

is a directed path in Hs, then so is 

(a(xi), . . . , a(x m )) — y (a(x 2 ), a(x m+1 )) — > ■■■ — > (a(x k+1 ), a{x k+m )) . 

We next show that there is an upper bound on the length of the directed path between any two 
vertices in the case when the graph Hs is strongly connected. 

Proposition 3.4. Suppose that Hs is strongly connected with period d. Then there is a positive 
integer N , a multiple of d, such that for any two points x and z in the same component Yi there is a 
path from x to z in the graph Hs of length N . Especially, between any two vertices in the graph Hs 
there is a directed path of length at most N + d— 1 . 

Proof. For two vertices x and y of the graph H, let D(x,y) denote the length of the shortest path 
from x to y. Observe that this is not a distance since it is not symmetrical in general. 

For each permutation ir £ & m pick a point x w in A^-, such that all its coordinates are greater 
than 1/2. Similarly, for all < i < d — 1 pick a point yi in Y^, such that all its coordinates are less 
than 1/2. Let K denote the maximum of the finite set 

{D(x n , yi) : A n CYi} U {D{y h x T ) : A T C Yi}. 

Note that K is a multiple of the period d. Let Qi denote the semigroup 

Qi = {k : there is a path from to of exactly length k} . 

Note that Qi C d ■ N. Furthermore, there is no multiple e of the period d such that Qi C e • N. Hence 
there exists a positive integer rrii such that d ■ N + mi C Qi. That is, there is path from yi to yi of 
any length which is a multiple of d and greater than or equal to mj. Let M be the maximum of mo 
through ma-i- Note again that M is a multiple of d. 

Let N = 2 • K + M. Pick two permutations ir and r such that A^, A r C Yi. We claim that there 
is a path from the point x n to the point x T of length N. We find this path by picking three paths. 
First, choose a path p\ from x^ to yi of length at most K. Second, choose a path p^ from to x T of 
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length at most K. Finally, pick a path p 2 from yi to of length M + 2 • K — D(x 7r , yi) — D(y^,x T ) 
which is at least M. By concatenating the three paths pi, p 2 and P3, the result follows. 

Now let show that there is a path from any point x in Yi to any other point z in Yi of length N. 
Let 7r and r be the permutations ir = H(x) and r = n(z). We do so using that we did choose 
the coordinates of x n and x T greater than the coordinates of y^. Namely, let e > be the smallest 
coordinate of the two points x and z. Let y be the point £ ■ yi- Now there is a monotone function a 
such that a(x 7T ) = x and a(yi) = y. Similarly, there is a monotone function j3 such that = y 

and /3(x T ) = z. Concatenating the three paths a(pi), e -pi and /3(ps) we obtain a path from x to z of 
length iV. □ 

Proof of TheoremUJh It follows from Proposition that T d : L 2 (Yi) — > L 2 (Yi). We will denote 
by A the restriction of T d to L 2 (lo)- Choosing a positive integer p with p ■ d > m we see that 
A p is a compact operator from L 2 (Yq) to itself. The operator A p has discrete spectrum which may 
accumulate only at so the same is true for A. We claim that, if tp is an eigenfunction of A with 
nonzero eigenvalue X d , then T*<£> is a nonzero eigenfunction of T d with eigenvalue X d and support in Y^. 
The only nontrivial part of this claim is that T l ip ^ 0. To see this, note that T d ~ l (T l ip) = X d ■ ip so 
T*(/? cannot be zero since A 7^ and tp is a nontrivial eigenfunction. 

Suppose now that A d G cr(A), let 93 be an eigenfunction of ^4 corresponding to the eigenvalue A , 
let uj be a dth root of unity, and let 



The function ip is nonzero because the right-hand terms are nonzero and have disjoint supports. A 
direct computation shows that Tip = oj ■ A • tp so the spectrum of T contains all of the numbers A • ui 
where X d G o~{A). On the other hand, any eigenvalue /i of T gives rise to an eigenvalue \i d of T d , so 
the nonzero spectrum of T consists exactly of the numbers oj • A where a; is a dth. root of unity and X d 
is an eigenvalue of T d . 

By combining Lemma 13.11 and Proposition 13.41 we have that for a non-negative, but non-zero, 
function / in L 2 (Y ) that T N (f) = A N l d {f) is a positive function. Hence the operator A is positivity 
improving and by Krem and Rutman |15t Theorem 6.3] the operator A has a positive spectral radius 
r(A). Furthermore, the spectral radius r(A) is a simple eigenvalue, all other eigenvalues are smaller 
in modulus and the eigenfunction p corresponding to r{A) is positive. Hence the spectral radius of T, 
r(T) = ^/r(A), is positive and a simple eigenvalue of T. Finally, the positive function tp is also an 
eigenfunction of T corresponding to the eigenvalue r(T). □ 

3.2 The directed graph Gs 

We examine relations between the infinite graph Hs and the finite graph Gs- We need the following 
lemmas, 

Lemma 3.5. Let x = (x±, . . . ,x m ) be a vertex in Hs such that U(xi, . . . ,x m ) = ir. Furthermore 
assume there is an edge in Gs labeled r leaving the vertex it. Then there exists x m+ \ in (0, 1) such 
that n(xi, . . . ,x m +i) = r. That is, there is an edge in Hs leaving the vertex x. 

Proof. Observe that r(m+l) is bigger than exactly r(m + l) — 1 of the numbers r(l), . . . , r(m). Hence 
pick x m +i such that it is bigger than exactly r(m + 1) — 1 of the numbers x\, — , x m . 

Iterating this lemma we obtain: 
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Lemma 3.6. Given a directed path from ir to a in the graph Gs- Let x be a vertex in Hg such that 
U(x) = tt. Then this directed path can be lifted to a directed path in Hs that ends with a vertex y such 
that II(y) = a. 

Now we give a sufficient condition for Hs to be strongly connected in terms of the graph Gs ■ 

Proposition 3.7. Let S C @ m+ i ; suppose that Gs is strongly connected, and suppose that the two 
monotone permutations 12 • • • (m + 1) and (m + 1) • • • 21 do not belong to the set S. Then the graph Hs 
is ergodic. 

Proof. We first prove that Hs is strongly connected. Let x = [x\, . . . , x m ) and y = (yi, . . . , y m ) be two 
vertices of Hs, and let tt = H(x) and a = H(y). Since Gs is strongly connected we can find a directed 
path from n to m ■ ■ ■ 21 in Gs- This directed path lifts to a directed path from x to z = (z±, . . . , z rn ) 
in where z\> z 2 > • • • > z m . 

For any directed graph G define the reverse graph G* to the graph G where we reverse the 
direction of each edge. In the reverse graph G* s we have a directed path from a to 12 • • • m. Lifting 
this directed path to a directed path in H* s and then reversing the path, we obtain a directed path 
from w = [w\, . . . , w m ) to y in Hs, where W\ < ■ ■ ■ < w m . 

Finally, there is a directed path from m ■ ■ ■ 21 to 12 • • • m in Gs- Hence we know that there is a 
directed path from u = (ui, . . . , u m ) to v = (v±, . . . , v m ) in Hs, where u± > • • • > u m and v\ < • • • < v m . 

Choose a so that < a < mm(z m , wi). We then have a directed path from a ■ u to a ■ v. Now, 
there is a directed path from z to a • u of length m, namely 

z = Oi, ...,z m ) — > (z 2 , ■ ■ .,z m ,aui) — > ■■■ — > (z m ,aui, . . . ,an m _i) — > (aui, . . .,au m ) = a-u 

using the fact that (m + 1) ■ ■ ■ 21 is not forbidden. We can now concatenate these five directed paths 
to obtain a path from x via z, via a • u, via a • v, via w, to y. 

Since Gs is strongly connected and has ml vertices, an upper bound on the length of the directed 
paths in Gs chosen above is ml — 1. Hence, the path that we have constructed has length at most 
3(m! - l) + 2m. 

To observe that Hs has period 1 note that we can construct a directed path from the vertex x 
to the vertex y that has length one more than the above construction. Namely, the path from tt to 
m ■ ■ ■ 21 can be extended by adding the loop (m + 1) • • • 21 at the end. Now by concatenating these 
two paths with a path from y to x we obtain two cycles whose lengths differ by one. Since the greatest 
common divisor of two consecutive integers is one, the graph Hs is ergodic. □ 

Recall that a permutation tt in 6 n is indecomposable if there is no index i such that 1 < i < n — 1 
and 7r(l), . . . ,ir(i) < i. Otherwise the permutation is decomposable. 

Proposition 3.8. Let S C (3 m+ i such that each permutation r G S is indecomposable. Then the 
graph Gs is strongly connected. 

Proof. Given two vertices tt = (tt\, . . . , TT m ) and a = (a±, . . . , a m ) of Gs- Since the integers tt\, . . . , TT m , 
a\ + m, . . . , a m + m are distinct, the following path (described by its edges) is well defined: 

n(7ri, . . . ,TT m ,o-\ + m),n(7r 2 , . . . ,TT m ,ai + m,a 2 + m), . . . ,n(7T m ,cJi + m, . . . ,a m + m). 

This path goes from the vertex tt to vertex a. Note that every edge r' on the path is decomposable. 
Hence this path avoids the forbidden indecomposable edges of S. □ 
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Proof of Theorem \1.5[ This follows directly from Theorems 11.21 and 11.71 and Propositions 13.71 and 13.81 

□ 

Proof of Corollary \l-4\ Let r be the single permutation in the set S. If t is one of the two monotone 
permutations the result will follow from descent pattern avoidance, see Theorem 14.21 If r is indecom- 
posable the result follows from Theorem 11.51 Finally if r is decomposable apply Theorem 11.51 to the 
upside down permutation {m + 2 — r(l), . . . , m + 2 — r(m + 1)) which is not decomposable. □ 

Note that there are examples of patterns S so that H$ does not have cycles even though the 
graph Gs has cycles. 

Example 3.9. S = {312,321}. In this case the graph Gs has a cycle. However, the graph H$ 
does not have a cycle and hence is not strongly connected. We observe this by noting that for a 
directed edge (x,y) — > (y,z) in Hs we have that x < max(y,z). Hence none of Xj's in a /c-cycle 
(x\, X2) — > (x2,X3) — > ■ ■ ■ — > (xk, x\) — > (xi,X2) can be the largest. 

Via the classical bijection tt 1 — > tt (see |19^ Section 1.3]) one obtains that the number {312,321}- 
avoiding permutations are in bijection with involutions, that is, permutations ir such that tt 2 = id. 
This was first observed by Claesson [3J. It follows that the generating function is exp(z + z 2 /2) and 
the asymptotic is 1/V2- exp(— 1/4) • (n/e) n / 2 • exp(y / n). 

On the other hand, if Gs does not have a cycle, then neither does H$- 

Lemma 3.10. Let S C © m +i and suppose that Gs has a directed cycle that contains the two vertices 
12 • • • m and m---21. Moreover, assume that the two monotone permutations 12- • • (m + 1) and 
{m + 1) • • • 21 do not belong to the set S. Then the graph Hs contains a directed cycle. 

Proof. Pick a vertex x = (x±, . . . , x m ) such that H(x) = 12 • • • m. The directed path from the vertex 
12 • • • m to the vertex m ■ ■ ■ 21 can be lifted to a path in Hs from the vertex x to a vertex y = 
(yi, . . . , y m ) where Tl(y) = m ■ ■ ■ 21. Similarly, pick a vertex z = (z±, . . . , z m ) such that H(z) = m ■ ■ ■ 21. 
The directed path from the vertex m ■ ■ ■ 21 to the vertex 12 • • • m can be lifted to a path in Hs from 
the vertex z to a vertex w = (w\, . . . , w m ) where IL(w) = 12- • • m. 

Using the two monotone functions a,/3 : (0, 1) — > (0, 1) defined by a(x) = (x + l)/2 and [3(x) = 
x/2, we have two directed paths: one from a(x) to a(y) and one from (3(z) to /3(w). 

Using that (m + 1) • • • 21 is an edge in Gs we have the following path from a(y) to 0(z), namely 

a (y) = {a(yi),---,a(y m )) — > (a(y 2 ), • . . , a(y m ), f3(zi)) — > > (P{zi), ■ ■ ■ , P{z m )) =/?(»• 

Similarly, we have the directed path 

P(w) = (f3(wi), . . . ,P(w m )) — > (f3(w 2 ), ■ ■ ■ ,(3(w m ),a(xi)) — > ■■■ — > (a(xi), . . . ,a(x m )) = a(x). 
Concatenate these four directed paths to obtain a directed cycle in Hs- □ 

4 Computational techniques 

In this section we discuss descent pattern avoidance which is a special case of pattern avoidance. First 
we introduce an analogue of the de Bruijn graph Djj, which has the advantage that it is smaller than 
the graph Gs- Moreover, if the graph Djj is ergodic so is the graph H S (jj) and we obtain that the 
associated operator is positivity improving. Second, for descent pattern avoidance we obtain that the 
eigenfunctions has a simplified form. Finally, in the last subsection we consider pattern avoidance 
that has symmetry. In these cases we show that we can obtain the adjoint eigenfunctions from the 
eigenfunctions. 
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4.1 Descent pattern avoidance 

The descent set of a permutation tt in the symmetric group on n elements is the subset of {1, . . . , n— 1}, 
given by {i : 7Tj > 7Tj+i}. An equivalent notion is the descent word, defined as follows. The descent 
word of the permutation tt is the word u(ir) = u\ ■ ■ -u n -i where U{ = a if tti < 7Tj + i and Uj = b 
otherwise. 

Let U be a collection of aft-words of length m. The permutation tt avoids the set U if there is no 
consecutive subword of the descent word of tt contained in the collection U. 

Descent pattern avoidance is a special case of consecutive pattern avoidance. For instance, per- 
mutations avoiding the word aab is the permutations avoiding the set S = {1243, 1342,2341}, since 
these three permutations are the permutations with descent word aab. More formally, for U a subset 
of {a, b} m define S(U) C & m+1 by 

S(U) = {tt £ & m+1 : u(tt)£U}. 

Then the set of permutations avoiding the descent words in U is the set of permutations avoiding 
S(U). 

For U a subset of {a, b} m define the associated de Bruijn graph Djj by letting the vertex set be 
{a,b} m ~ l . For x,y £ {a, 6} and u E {a, b} m ~ 2 such that xuy U let there be a directed edge from 
xu to uy. When the set U is empty, the graph D\j is the classical de Bruijn graph D m ^\. 

Lemma 4.1. Let U be a subset of {a, b} m . If there is a cycle c of length N > n + 1 that do not consists 
only of the loop a m or not only of the loop b m , then the cycle c can be lifted to a cycle of length N in 
the graph Hgnjy 

Proof. Let c be the cycle 

VlV2 ■ ■ ■ V m -1 ► -"2^3 ■■■V m >-" > V N Vi ■ ■ ■ V m - 2 ► ^2 • • • 

where each Vi is either a or b. We would like to pick N real numbers x±, . . . ,xn in the interval (0, 1) 
such that 

Xi < x i+ i if Vi = a, , . 

Xi > x i+ i if vi = b, 

where all the indices are modulo N. Since all the letters v\ through vn are not the same, we may 
assume that vn-i = a and vn = b. Pick x\ arbitrarily. Pick x 2 through xn-i such that inequality (|4.ip 
is satisfied. Finally, pick x^ in the interval (max(xjv-i,^i)) 1)- Now in the graph Hgnj\ we have the 
cycle 

(xi,X 2 , Xm) > (X2,X 3 , . . . , Xm+l) > ' ' ' > (XN, %l, ■ ■ ■ , X m -l) > (xi,X 2 , ■■■ , X m ). 

□ 

Theorem 4.2. Let U be a subset of {a,b} m . If the de Bruijn graph Djj is strongly connected, then 
the graph Hs(U) * s a ^ so strongly connected. Furthermore, the de Bruijn graph Djj has the same period 
as the graph Hgnn . 

Proof. Note that the graph Djj has the directed edge a m ~ 1 b since otherwise it would not be strongly 
connected. Given two vertices x and y in Hgmy To prove that Hgnj\ is strongly connected it is 
enough to find a directed path from x to y. 
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We can find a path from u(U(x)) to a m_1 in the graph Djj that consists of at least m + 1 edges. 
Similarly to the lifting lemma, Lemma |3.6| we can lift this path to a path in Hs(jj) that starts at the 
vertex x and ends, say, in the vertex z = (z±, . . . , z m ). Note that Z\ < ■ ■ ■ < z m . Moreover, we can 
find a path from a m ~ 2 b to u(H(y)) in Djj that has length at least m + 1. Lift this path to a path that 
ends in the vertex y and begins at w = (w\, . . . , w m ), where w\ < ■ ■ ■ < w m ^\ > w m . 

Let Vi = max(zj, for 2 < i < m. Observe that we have the string of inequalities z\ < V2 < 

■ ■ ■ < v m > w m . We can now concatenate these two paths as follows. Replace each occurrence of Z{ 
and Wi-i by Vi for 2 < i < m in each of the two paths. Then we may connect the vertex (zi,V2, ■ ■ ■ , v m ) 
with the vertex (v2, ■ ■ ■ ,v m , w m ) via the edge that goes across the edge with descent word a m ~ 1 b. Thus 
the graph H S (u) is strongly connected. 

Since there is a graph homomorphism from Hg^jj-j to D\j we know that the period of D\j divides 
the period of H S (jj)- To see that the periods are equal, pick a vertex w of Djj that differs from a m ~ x 
and b m ~ 1 . Then any cycle of length greater than m + 1 through the vertex w in Djj lifts to a cycle 
of the same length in Hg^jy Hence the greatest common divisor of lengths of cycles through w is 
a multiple of the greatest common divisor of the cycle lengths of Hgmy Hence the two periods are 
equal. 

Note that this argument only works when m > 3 since there is no such vertex w in the m = 2 
case. But the remaining m = 2 case is straightforward to check. □ 

4.2 Invariant subspace for descent pattern avoidance 

For an a6-word u of length m — 1 define the descent polytope P u to be the subset of the unit cube [0, l] m 
corresponding to all vectors with descent word u. That is, 

P u = {(xi,. . . ,x m ) £ [0, l] m : Xi < x i+ i if Ui = a and x { > x i+ \ if m = b}. 

Observe that the m-dimensional unit cube is the union of the 2 m_1 descent polytopes P u . Now the 
operator T corresponding to the descent pattern avoidance of the set U has the following form. For 
an a6-word u of length m — 2 and y £ {a, b} we have 

T(f)\ Puv = H X(auy)- f(t,x 1 ,...,x m „ 1 )\ PaK dt (4.2) 
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+ / X{buy) ■ f(t,xi,...,x m -i)\ P dt, 



where by abuse of notation we let x( w ) = 1 if w does not belong to the set U and x{ w ) = otherwise. 

Proposition 4.3. Let T be the operator associated with a descent pattern avoidance and k is an 
integer such that < k < m — 1. Let u be an ab-word of length m — 1. Then the function T k (f) 
restricted to the descent polytope P u only depends on the variables x\ through x m _/ t . 

Proof. Proof by induction on k. When k = there is nothing to prove. When 1 < k < m — 1, we know 
by induction that the restriction of T k ~ 1 (f) only depends on x±, . . . , x m _k+i- By the shift of variables 
in the right hand side of equation (|4.2p . we obtain that T k (f) does not depend on the variable x m -k+i, 
completing the induction. □ 

Corollary 4.4. Let T be the operator associated with a descent pattern avoidance and let ip be an 
eigenfunction associated with a non-zero eigenvalue A. Then the eigenfunction restricted to each 
descent polytope P u only depends on the variable x\. 
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Proof. Since X m ■ tp = T m (tp) the eigenfunction has the required form. □ 

Corollary 4.5. Let T be the operator associated with a descent pattern avoidance and let tp be an 
eigenfunction associated with a non-zero eigenvalue A. Assume that f is a generalized eigenfunction, 
that is, it satisfies the equation A • / = T(f) + tp. Then the function f restricted to each descent 
polytope P u only depends on the variable x\ . 

Proof. By induction on k. Assume that 1 < k < m and that / restricted to each descent polytope only 
depends on the variables x\ through x^. Then T(f) + tp only depends on x\ through x^-i showing 
that A • / only depends on x\, . . . , xt-i- □ 

Let V be the subspace of L 2 ([0, l] m ) consisting of all functions / that only depend on the variable x\ 
when restricted to each of the descent polytopes P u . Observe that the subspace V is invariant under 
the operator T. That is, the operator T restricts to the subspace V. Moreover the constant function 1 
belongs to V. Hence to understand the behavior of T n (l) it is enough to study this restricted operator. 

In order to describe the subspace V more explicitly define for an aft-word u of length m—1 the 
polynomial h{u]X\) as follows: 



h(u;xi)= / ldx2 ■ ■ ■ dx m . 

J (xi,X2,..;X m )eP u 

These polynomials were first introduced and studied in [8], with different notation. 

Let p be a vector (p u (%i)) ue { a M.m-i- That is, the vector p consists of one-variable functions in the 
variable x\ and is indexed by ah- words of length m—1. Consider the function / on [0, l] m defined by 

f{x 1 ,...,X m )\p u =p u (x 1 ) 

for all afc-words u of length m—1. Observe that the function / belongs to L 2 ([0, l] m ), and hence to 
the invariant subspace V, if and only if 

l 

h(u;xi) ■ \p u (xi)\ 2 dx\ < oo 



o 



for all ab- words u of length m — 1. For two functions / and g in the subspace V, corresponding to the 
two vectors (p u (xi)) ue ^ b y m -i and (q u (xi)) ue ^ b y m -i, the inner product is given by 



(fi9)= "52 / h(u;xi) -p u (xi) ■ q u (x\)dxi. 



We end this section by a structural result about the subspace V. 

Proposition 4.6. The invariant subspace V is isometrically isomorphic to the Hilbert space 

Proof. The isomorphism of the Hilbert spaces V — > L ([0, 1]) is given by 

(Pu(a?l))« e {o,6}»»-i >• (VHuiX!) -pu(xi)) . 

□ 
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4.3 Symmetries 

Let J and R be the following two involutions on the space L 2 ([0, l] m ): 

(J/)(xi,x 2 , • • • ,x m ) = f(l - x m , . . . , 1 - x 2 , 1 - Xl), 

(Rf){x\,X2, . . . ,X m ) = f(x m ,...,X 2 ,Xi). 

Observe that both J and R are self adjoint operators. 

Lemma 4.7. Assume that x has the symmetry 

x(x\,x 2 , ■ ■ ■ ,x m ,x m+ i) = x(l - x m+ i, 1 - x m , . . . , 1 - x 2 , 1 - xi). 

Then the adjoint of the associated operator T is given byT* = JTJ. Moreover, ifip is an eigenfunction 
of the operator T with eigenvalue A then ip = Jip is an eigenfunction of the adjoint T* with the 
eigenvalue X. Furthermore, we have the equality (l,V0 = (v 9 ;!)- 

Proof. We have that 

JTJf(xi,x 2 ,...,x m ) = JT/(1 - x m , ■ ■ ■ , 1 - x 2 , 1 - xi) 



J / x(t, xi,..., x m ) ■ f(l - x m -i, . . . , 1 - xi, 1 - t)dt 
Jo 

/ x(t, 1 - x m , . . . ,1 - xi) ■ f(x 2 , x m , 1 - t)dt 
Jo 

/ X(l - t, 1 - x m , . . . , 1 - xi) ■ f(x 2 , ...,x m , t)dt 
Jo 



X(xi, ...,x m ,t)- f(x 2 , x m ,t)dt 
T ■ ■ • j x m — i, x m ). 







For the second statement consider the following line of equalities T* Jip = JTJ Jip = JTip = A • Jip. 
Lastly, = (1,7$) = (l,Jlp) = {Jl,Tp) = {l,Tp) = {ip,l). □ 

Similarly to Lemma 14.71 we have the next lemma. Its proof is similar to the previous proof and 
hence omitted. 

Lemma 4.8. Assume that x has the symmetry 

X(XI,X 2 , . . .,X m ,X m+ i) = x(x m +l,X m , ■ ■ ■ ,X 2 ,Xl). 

Then we have that the adjoint of the associated operator T is given by T* = RTR. Moreover, ifip is an 
eigenfunction of the operator T with eigenvalue A then tp = Rip is an eigenfunction of the adjoint T* 
with the eigenvalue A. Furthermore, we have the equality = (v 9 )!)- 

5 123- Avoiding permutations 

A 123-avoiding permutation is a permutation ir G & n with no index j so that ttj < vrj + i < TTj+ 2 , where 
1 < J ' < n — 2. Let a n (123) denote the number of 123-avoiding permutations in <5 n . 
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5.1 Eigenvalues and eigenf unctions 

Since 123-avoiding permutations can be viewed as permutations with no double descents Corollary 14.41 
allows us to recast the problem of finding eigenfunctions in two variables into finding two one-variable 
functions. 

Proposition 5.1. The eigenvalues of the operator T are given by 



V3_ 

2-7T-(fc + i) 



X k = - _ — , lv (5.1) 



where k G Z and the associated eigenfunctions are given by 

x cos(f + f -f) <fO<x<y<l, 

tp k = exp -— - ■ { ) ( (5.2) 



2 • A 7 I sin ( f + ^ • § ) £/ < y < x < 1. 



Proof. Avoiding the pattern 123 is equivalent to avoiding the descent set pattern aa. Hence Corol- 
lary H3J states that the eigenfunctions <p can be written as 



p(x) if < x < y < 1, 
g(jc) ifO<y<x<l. 



Then the defining equations for eigenvalues and eigenfunctions reduces to the integral system: 

X-p{x) = [ q(t)dt, (5.3) 



/ p(t)dt+ [ q(t)dt. (5.4) 
Jo Jx 



A • q(x 

First, differentiating with respect to x, we obtain the first-order system 

X-p(x) = —q(x), (5-5) 
X-q'(x) = p(x) — q(x). (5-6) 

These equations have only the trivial solution if A = 0, so A = is not an eigenvalue. If A ^ then 
the first-order system (I5.5p ~ (l5.6p implies the second-order equation 

A 2 • p"(x) + A • p'{x) + p(x) = 0. 

This equation has the general solution 

oj \ ( 



p(x) = A ■ exp ( — • xj + B ■ exp ( — • x ) , (5.7) 



where u = exp (^-f-^)- That is, ui satisfies the relation oo 2 + to + 1 = 0. Moreover, equation ()5.5 
implies that 



-to ■ A - exp ^— • x*J — to 2 • B ■ exp ( — ■ x ) . (5. 
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Setting x = and x = 1 in equations (|5.3|) and (|5.4|) and using that A/Owe obtain the boundary 
conditions: 

p(0) = g (0), (5.9) 
p(l) = 0. (5.10) 

Substituting the expressions for p(x) and q{x) from equations (|5.5p and (|5.6p into boundary condi- 
tion (|5.9|) we obtain A + S = — cj • A — u 2 ■ B. This is equivalent to u ■ A + B = 0. Hence we may 
set A = 1/2 • exp (^) and £? = A = 1/2 • exp (— Substituting equation ()5.5|) into the second 
boundary condition (|5.10p implies that 



A ■ exp ^— ^ = —5 • exp = u • A - exp 



Cancelling A on both sides and taking the logarithm gives 

UJ 2 • 7T • i UJ 1 

where is an integer. Since cj — w 2 = v3 • i we obtain expression (|5.ip . Moreover p(x) is given by 

1 (n;-i\ ( x\ 1 / 7T-A / 2 x 

= --exp^— j-exp^- J + --exp^-— j-exp^ ■- 

/ X \ / 7T \/3 X \ 

Now equation (|5.5p implies the claimed expression for g(x). □ 

Note that the eigenvalues are ordered by 

A > -A_i > Ai > -A_2 > A 2 > -A_3 > A 3 > • • • > 0. 

Furthermore, the calculations in Proposition 15. II showed that they all have a unique eigenfunction. It 
remains to show that they are simple. 

Proposition 5.2. The eigenvalues A of the operator T are simple, that is, they have index 1. In other 
words, there is no function f(x,y) such that \k ■ f = T(f) + (pi-. 



r(x) if < x < y < 1, 
s(x) if < y < x < 1. 



Sketch of proof. Using Corollary 14.51 we can write the function / as 

f(x,y) = 

Then equation reduces to the integral system: 

A • r{x) = [ s(t)dt + p(x), 
A • s(x) = 



r(t)dt+ I s(t)dt + q(x). 

J X 
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Note that we obtain the two boundary conditions r(l) = and r(0) = s(0). Next differentiating with 
respect to x, we have the first-order system: 

A • r (x) = —s(x)+p'(x), 
X-s'(x) = r(x) — s(x) + q'(x). 

This system of differential equations can be solved as in the proof of Proposition 15.11 However, the 
solution does not satisfy the boundary conditions, completing the sketch. 
By applying the involution J we obtain the adjoint eigenfunction 

cosf^ + ^-i^O if 0<x <y<l, 



2-A / sinf^ + ^-i^) ifO<y<x<l 



Proposition 5.3. Let A be an eigenvalue of T with eigenf unctions (p and let ip be the eigenfunction 
of the adjoint operator T* with eigenvalue A. Then the following identities hold: 

(¥>,!) = (1,?) = ^-A 2 , (5.12) 

faty = |.(-l)*.A.«p(-iV (5.13) 

In particular 



(*>,?) 1 ' ' V2-A 

Proof. In the following calculations we use the facts that cos(\/3/(2 • A)) = (— l) fc /2 and sin(\/3/(2 • 
A)) = (— l) fc \/3/2. We also use the expression for ip in equation (|5.2p . First, we note that 

(cp, 1) = / p(x)dxdy + / q(x)dxdy 

J0<x<y<l J0<y<x<l 
■1 H 

(1 — x)p(x) ■ dx + x • q(x)dx. 
Jo Jo 

Explicit computation shows that 

(l-x)-p(x)dx = ^.A 2 -(l-(-l) fc .exp(-l/(2-A))) 

x-q{x)dx = ^-A 2 -(-l) fc -exp(-l/(2-A)) 
which shows (|5.12[) . Next, using (|5.2p and (|5.1ip and dropping subscripts as before, we have 



(<p,ip) = / p{x) ■ p(l - y)dxdy + q(x) ■ q(l - y)dxdy 

J0<x<y<\ J0<y<x<l 

= J (p(x) - J p(i - y)dy + q{x) - J <?(i - y)dy^j dx. 

Carrying out the y integration and simplifying, we obtain 

{<p,V) = ^-(-l) k - j\-exp (-^rj) dx 
which gives f)5. 13|) . □ 
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5.2 Asymptotics 

The above computations show that all eigenvalues of T are simple and give the eigenvalues and the 
coefficients explicitly. We thus obtain the following expansion for (T n_2 (l),l) = a n (123)/n! as an 
immediate consequence of Theorem 11.14 Propositions 15,11 and 15.31 

Theorem 5.4. Let K be a non-negative integer. The number of 123-avoiding permutations satisfies 
the following asymptotic expansion 



corresponding to the 2K + 1 largest eigenvalues. 

6 2 13- Avoiding permutations 

A 213-avoiding permutation is a permutation tt £ & n which contains no sequence of the form 

TTj +1 < TTj < TTj +2 

for any j with 1 < j < n — 2. We denote the number of 213-avoiding permutations of <3 n by a n (213). 
Thus, S consists of the single permutation 213 and 



By symmetry, the study of 213-avoiding permutations is equivalent to 132-avoiding permutations, 
231-avoiding permutations and 312-avoiding permutations. However the case of 213-avoiding permu- 
tations gives the most straightforward equations. 

6.1 Eigenvalues and eigenfunctions 

In what follows, we will make use of the error function 




where is given 




sum contains 2K + 1 terms 





(6.1) 



which extends to an entire function on C, and the function 




(6.2) 



Let 





Now we characterize the nonzero eigenvalues and eigenfunctions. 
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Proposition 6.1. The non-zero eigenvalues A of the operator T satisfy the equation 

and the corresponding eigenf unctions are 

' Q( x )~J- nq(t) dt ifO<x<y<l, 
q(x) if < y < x < 1, 



where q(x) is given by A6.2\) . 

Proof. The defining relations for the eigenfunctions are 



X-p{x,y) = [ p(t,x)dt+ [ q(t,x)dt, (6.4) 

JO Jy 
rx rl 

X-q(x,y) = / p(t,x)dt+ / q(t,x)dt. (6.5) 

JO Jx 

Now observe that in the right-hand side of equation (|6,5p there is no dependency on the variable y. 
Hence we may replace q(x,y) with q(x). Now subtract equation (]6.5p from equation (|6.4|) 

rv 

A • (p(x, y) - q(x)) = - I q(t)dt. 



That is, 

1 [ y 

p[x,y) = q(x) - - J q(t)dt. (6.6) 



Substitute equation (|6.6p into equation (|6.5D : 



A • q(x) = J (q(t) - i J q{s)ds^j dt + j q(t)dt 



q(t)dt — - • / / q(s)dsdt 
o A J J t 

P 1 PX f k S 

/ q(t)dt - - ■ / / q(s)dtds 
Jo * J J 



q(t)dt — — ■ j s ■ q(s)ds. 



Hence we have the following integral equation for q{x) 

rl rx 

A 2 • q(x) = A • / q(t)dt - s- q(s)ds. (6.7) 
Jo Jo 

Differentiating once we have 

A 2 • q'(x) = —x ■ q(x). (6.8) 
The solution to this differential equation is 

q( x ) = C-e X p(~^\. (6.9) 
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By setting the constant C to be 1 we obtain a solution to equation (|6.8|) . Now substitute this solution 
for q(x) into the integral equation (|6,7p and set x = 0: 



A 2 = A • / exp ) '// 



o 



2 • A 2 

l/(\/2-A) 



V2 • A • / exp (-ti 2 ) 

v^-A 2 



^_ erf(l/(V2-A)), 

where the substitution in the integral is u = t/{y2- A). Hence the non-zero eigenvalues A satisfy 
equation ()6.3|) . □ 

Proposition 6.2. The eigenvalues A of the operator T that have index greater than or equal to 2, 
satisfy the equation 

exp ( -—"-2 J = A - 1. (6.10) 



2- X\ 

In other words, if there is a function f(x,y) such that A • / = T(f) + tp this implies that A satisfies 
equation \6.10}) . 

Proof. The defining relations for the eigenfunctions are 

px pi 

\-r(x,y) = / r(t,x)dt+ / s(t, x)dt + p(x,y), (6-11) 
Jo Jy 

X-s(x,y) = / r(t,x)dt+ / s(t, x)dt + q{x). (6.12) 

JO Jx 

Note that in the right-hand side of equation (16. 12H there is no dependency on the variable y. Hence 
s(x,y) is a function of x only, and we write s(x) henceforth. Subtracting equation ()6. 12j) from equa- 
tion (I6.1ip and dividing by A gives 

I rv i 
r(x,y) = S ( X )-J'J s (t)dt + (p(x,y) - q{x)) 

= s(x) - ^ • j\\ ■ s(t) + q(t))dt, 

where the last step is by equation ()6.6p . Substituting this expression into equation ()6. 12j) yields the 
integral equation for s(x): 



1 



A • s(x) = I ( s(t) - — • / (A • s(u) + q{u))du )dt+ / s(t)dt + q(x) 



1 

A 2 " 
1 

a 2 " 
1 

A 2 
1 

A ' 



(A • s(u) + q(u))dudt + / s(t)dt + q(x) 
(I Jt Jo 

x pu pi 

/ (A • s(u) + q(u))dtdu + s(t)dt + q(x) 
(I Jo Jo 

1 



u ■ (A • s(u) + q{u))du + I s(t)dt + q(x) 

Jo 

[ u-s(u)du+[ s(t)dt-\- f q(t)dt + 2-q(x), (6.13) 
Jo Jo A 7o 
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where the last step is using equation (|6.7|) . Differentiate with respect to x to obtain the first order 
differential equation in s(x): 

A • s (x) = — — ■ x ■ s(x) + 2 • q'(x). 
A 

The general solution to this differential equation is 

x^ ( x^ \ 

4 X ) = ~J2 " eX P ( ~2~\2 J + C ' 

where we used that q{x) = exp ^— -^jj- Observe that the homogeneous part of this solution is 

expected. It corresponds to the homogeneous part of the equation A • / = T(f) + (p. Hence we may 
set C = without loss of generality. Setting x = in (|6.13p yields 



0=/ s{t)dt-\- [ q(t)dt + 2- q(0). 
Jo A Jq 

Using that q(0) = 1, Jq 1 q(t)dt = A and 



1 \ h - „ / 1 



s(t)dt = exp — - — — — a — ■ A • erf 



2-AV V 2 \V2-X 

we obtain the desired equation. □ 

For completeness we state: 
Lemma 6.3. A function in the kernel of the operator T has the form 



P(x,y) if0<x<y<l, 
if0<y<x<l, 



where p{x,y) satisfies f x p(t,x)dt = 0. 
The adjoint operator T* is given by 

Jo <?(y, u)du + jy p(y, u)du if < x < y < 1, 
Jo q(y, u ) du + Jy P{y, u)du if <y < x < 1. 



T*(f(x,y)) 



Proposition 6.4. For a non-zero eigenvalues A of the operator T the corresponding eigenfunction of 
the adjoint operator T* is 



i>{x,y) 

where 



P*(y) ifi<x<y<l, 

p*(y) - j ■ IxP*( u ) du ^ <y < x < 1, 

,2 \ / „,2 



p*(y) = -2 ■ y ■ exp U^) + 2 • A + • y ■ exp U^) • erf . (6.14) 
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Proof. The defining relations for the eigenfunctions are 

ry 



X-p*(x,y) 
A • q*{x,y) 







*(y,u)du+ / p*(y,u)du, 



*(y,u)du+ / p*(y,u)du. 



(6.15) 
(6.16) 



Observe that there is no dependency on the variable x in equation (|6.15p . Thus we write p*{x,y) 
p*(y). Subtracting these two equations we have 



A- (q* (x,y) -p*(y)) = - / p*{u)du, 

J X 



such that 



f{x,y)=p*(y)---l p*{u)du. 



(6.17) 



Substituting this expression into equation (|6.15p one obtains 



X-p*{y) 



p*{u) — — ■ I p*(v)dv J du+ I p*(u)du 



p*{u)du — — 
A 



1 f y f 1 



Jy 



p* (v)dvdu 



p*(u)du— — -y- I p*(v)dv. 



That is, p*(y) satisfies the integral equation 

X 2 -p*(y) = X- / p*{u)du — y- I p*(v)dv. 
Jo Jy 

Differentiating this equation twice we obtain 

A 2 -y>*'(y) = y-p*(y)- p*(v)dv, 

Jy 

A 2 -P*"(y) = yp*'(y) + 2-p*(y). 
The solution of this differential equation is given by 

,2 



P* (y) = Ci ■ y ■ exp 



y 



2- A 2 



+ c 2 - 



2 • A + V2 • vr • y ■ 



cxp 



2 • A 2 



• erf 



a/2 • A 



(6.18) 

(6.19) 
(6.20) 

(6.21) 



Setting y = in equations (|6.18p and (|6.19p we obtain A • p*(0) = Jq p*(u)du = — A 2 • p'(0). Inserting 
this condition into the solution of the differential equation (16.21H we obtain C\ = — 2 ■ C%. Moreover 
setting Ci = 1 we obtain equation (|6.14p . □ 



Lemma 6.5. A junction in the kernel of the adjoint operator T* has the form 

ip(x,y) = 

where q*(x,y) satisfies Jq q*(y, u)du = 0. 



if0<x<y<l, 
Q*(x,y) if0<y<x<l 
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Proposition 6.6. For a non-zero eigenvalue A with eigenfunction ip and eigenfunction tjj of the adjoint 
operator T* , we have 

(¥>,!) = A 2 , 
(1,^) = 2 • A 3 , 

(<p,tf) = 2-A 2 -exp(-l/(2-A 2 )). 



In particular, 



{(p,l) ■ (1,V) , :i . ,, , 2 , 



A • exp(l/(2 • A )). 



Proof. In the calculations that follows we will use the relations erf(l/(\/2- A)) = \/2/tt, Jq q(x)dx = A 
and ^p*(y)dy = 2-X 2 . 

First the inner product between the eigenfunction and the constant function V. 

((p, 1) = J q(x)dxdy — — ■ j f q(t)dtdxdy 

J[0,1] 2 A Jo< x <y<lJ x 

= [ q(x)dx -\ - ( t ■ (1 - t) ■ q{t)dt 
Jo A J 

= A-A + ^-A 2 -erf(l/(V2-A)) 
v2 

= A 2 . 

Second, the inner product between the adjoint eigenfunction and the constant function 1. We have 
(l,Y>) = / P*{y)dxdy — — - j j p*(u)dudxdy 

J[0,1] 2 A Jo<y<x<\ Jx 

f 1 1 f 1 u 2 

= J p*(y)dy-j-J —-p*{u)du 



2^A ' ' ;! [jl)ihi 



-ro 



The first integral is given by ii = (A - 2 • A 2 - 2 • A 3 ) • exp(l/(2 • A 2 )) - 1/3 + 2 • (A + A 2 + A 3 ). The 
second integral we solve by integration by parts letting /' = ^1 — • \/2 • tt • y • exp anc ^ 

9 = erf (tIa) • Then So ^ gdy = ^0 ~ So f 9 ' dy is given by: 

r.Af „ „ ^ I 1 

h = 



j= ■ A • (2 • A + 2 • A 2 - y 2 ) ■ exp(y 2 /(2 • A 2 )) • erf (y/(V2 • A)) 



J 



" /< 2 

Jo 



A + 2 • A 2 - y 2 )dy. 
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Combining all the terms in the sum I\ + I2 we obtain 2 • A 3 . The third inner product is given by 



(<P,W) = / q{x) ■ p* (y)dxdy 
J[o,i} 2 

— — • f f q(t)dt ■ p* (y)dxdy — — ■ f q(x) ■ f p*(u)dudxdy 

A Jo<x<v<l Jx A Jo<v<x<l Jx 



IO<x<y<Ux * J0<y<x<l 

1 \ / ,1 x 2 



q(x)dx • / P*(y)dy - T • / t ■ q(t) ■ p*(y)dtdy 

J \Jo J A J0<t<y<l 

2 • A 3 - 2 • A • f (1 - exp(-y 2 /(2 • A 2 ))) • /(y)^ 
Jo 

-2 • A 3 + 2 • A • / exp(-y 2 /(2 • A 2 )) • 
Jo 

-2 • A 3 + 4 • A • ^ (-y + exp(-y 2 /(2 • A 2 )) + yft/ V2 ■ y ■ erf (y/( V2 ■ A))) dy 



2-A 2 -exp(-l/(2-A 2 )). 



□ 



6.2 Asymptotics 



We know that the largest root Ao of the eigenvalue equation (|6.3[) is real, positive and simple, since 
the associated operator T is positivity improving. However, to say a bit more about the eigenvalues 
consider the related equation eii(z) = \/2/tt. 



Since the error function is an increasing function on the real axis, the equation erf (z) = \/2/tt 
has a unique real root. The error function is an odd function hence we know by the strong version 
of the little Picard theorem that the equation erf(z) = y2/vr has infinitely many roots. Moreover, 
the complex roots appear in conjugate pairs. To summarize this discussion we have: The eigenvalue 
equation has a unique real root which is positive and is the largest root. The remaining infinitely 
many roots are all complex and appear in conjugate pairs. 

Numerically, we can approximate the roots of the eigenvalue equation (|6.3p . The unique real root 
is Ao = 0.7839769312 .... The next four largest roots are: 

Ai, 2 = 0.2141426360 . . . ± 0.2085807022 ... -i 
A34 = -0.1677323922 . . . ± 0.2418627350 ... -i 

Furthermore, these five roots do not satisfy equation (I6.10p in Proposition 16. 2^ so they are simple 
eigenvalues. From Proposition 16.61 we conclude: 

Theorem 6.7. The number of 213 -avoiding permutations satisfies 

a„(213) _ / 1 



6XP (2 J A 2 )" A ° +1 + ° (|Al|n) 



where Ao = 0.7839769312. . . is the unique real root of the equation erf (l/(\/2 • A)) = \f2/ir, and Ai 
is the next largest root and its modulus is given by \\\\ = 0.298936411 .... 

By considering the next two conjugate roots Ai and A2 we obtain an approximation for the next 
real term in the expansion of o n (213)/n!: 

2 • 1.158597034 . . . • (0.2989364111 . . .) n+1 • cos (-5.593221320 . . . + (n + 1) • 0.7722415374 . . .) . 
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7 123, 231, 312- Avoiding permutations 

Let us consider 123, 231, 312-avoiding permutations. Now the set S is not a singleton, but has cardi- 
nality three, that is, S is given by the set {123,231,312}. The function % is given by 



X(xi,x 2 ,x 3 ) 
That is, the operator T is described by 



1 if X3 < X2 < Xi, 

1 if X2 < X\ < X3, 

1 if X\ < X3 < X2, 

otherwise. 



(Tf)(x,y) 



%f(t,x)dt X0<x<y<l, 
jy f(t, x)dt + £ f(t, x)dt if < y < x < 1. 



7.1 Eigenvalues and eigenfunctions 

For a real number t let {t} denote the fractional part of the real number t, that is, {t} = t — [t\. 
Observe that the fractional part belongs to the interval [0, 1). 

Lemma 7.1. The subspace W of L 2 ([0, l] 2 ) consisting of functions of the form 

f(x,y)=g({x-y}) for (x, y) G [0, l] 2 (7.1) 
is invariant under the operator T. 

Proof. There are two cases to consider. First, if x < y, we have 

T(g({x - y})) = [ V g({t - x})dt = T g(s)ds. 

Jx JO 

Second, if y < x we have 

T(g({x - y})) = [ V g({t - x})dt + f g({t - x})dt 

JO Jx 
rl+y—x rl — x rl+y—x 

= / g(s)ds+ / g(s)ds = / g(s)ds. 

Jl-x Jo Jo 



These two cases can be summarized as the integral from to {y — x}. However the upper bound of 
the integral can be written as 1 — {x — y}, that is, the operator can be expressed as 

l-{x-y} 



ri-^x-y} 

T(g({x - y})) = / g(s)ds. 
Jo 



□ 



Since constant function 1 belongs to W, it is enough to consider eigenvalues and eigenfunctions in 
the subspace W. 
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Proposition 7.2. The eigenvalues of the operator T restricted to the invariant subspace W are given 
by Xj = (— 1)0' —1 )/ 2 ■ 2/(tt ■ j) where j is an odd positive integer and the corresponding eigenf unctions 
are described by 

cp(x,y) = cos({x - y}/\). 



Proof. The defining relations for the eigenfunctions are 



-l-i 

A • g(t) = I g(s)ds. (7.2) 
o 



Observe that A = is not an eigenvalue. Differentiating equation (j7.2|) twice, we obtain 

X-g'(t) = -g(l-t), (7.3) 
A •</'(*) = g'(l-t) = -l/X-g(t). (7.4) 



This second order equation has the general solution A ■ cos(i/A) + B ■ sin(t/A). Setting t = 1 in (|7.2p 
implies that g(l) = and then setting t = in (|7.3|) implies that g'(0) = 0. Hence B = and the 
eigenfunction is cos(t/A). For this function, equation (|T.2[) implies that cos(t/A) = sin((l — t)/X). 
Setting t = we obtain that the eigenvalues satisfy sin(l/A) = 1. They can be parametrized as in the 
proposition. □ 

Note that the eigenvalues tend to as 

Ai > -A 3 > A 5 > -A 7 > • • • . 

Proposition 7.3. The eigenvalues X of the operator T restricted to the subspace W are simple. That 
is, there is no function h(t) such that X ■ h = T(h) + g where g(t) = cos(i/A). 

Proof. The defining relations for the eigenfunctions are 

,i-t 

X-h{t) = / h(s)ds + g(t). (7.5) 



o 



Differentiating this equation twice gives 



X-h'(t) = -h(l-t) + g'(t), (7.6) 
A • h"(t) = h'{l-t)+g"{t). (7.7) 

Substituting (17.61) into (17. 7h and using (17. 4p . we obtain the second order differential equation: 

X 2 -h"{t) + h{t) = 2- X-g"(t). 

The general solution is 

hit) = -4r • t ■ sinit/X) + Ci • cos(t/A) + C 2 • sin(t/A). 
X z 

In the following we use that sin(l/A) = 1 and cos(l/A) = 0. Observe that setting t = 1 in ()7.5[) gives 
X-h(l) = g(l), which implies C 2 = 1/A 2 . Similarly, setting t = in ([75]) gives X-h'(0) = -h(l)+g'(0) 
which implies C2 = 1/(2 • A 2 ), a different value for C 2 . Hence there is no such a function h satisfying 
the integral equation (I7.5p . □ 
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Observe that the function x satisfies the symmetry in Lemma [4 .71 Hence the adjoint eigenfunction 
tp is given by J(p. However, since the value of ip(x, y) only depends on the difference x — y, we have 
that Jip = ip. 

Proposition 7.4. For a non-zero eigenvalue A with eigenfunction ip and eigenfunction ip = if of the 
adjoint operator T* , we have 

(</>,l) = (1,?)=A, 
(<p,ty = 1/2. 



In particular, 



(^i)-M) _ 2 A2 



7.2 Asymptotics 

The above computations show that all eigenvalues of the operator T are simple and given explicitly. 
We thus obtain the following asymptotic expansion for (T n ~ 2 (l),l) = a n (123, 231, 312)/n! as an 
immediate consequence of Theorem 11.11 Propositions 17.21 and 17.41 



" " ' ' 23 ' 231 ' 312) - = 2 • ]T Xj + O (|A 2fc+3 D • 



j odd 



However, when we let k tend to infinity, observe that the right hand side converges. Hence we obtain 
the exact expression 

q n (123,231,312) ^ u N-i)/2-n ( 2 Y 

— 5 — h 

j odd 

for n > 2. However this is the expression for E n _i/(n — 1)! where E n denotes the nth Euler number; 
see equation (|1 .8j) . Thus we conclude that 

Theorem 7.5. The number of 123, 231, 312 -avoiding permutations in & n is given by n ■ E n -% for 
n > 2. 

This result is due to Kitaev and Mansour [13J. Note that the special form of the eigenfunctions and 
the invariant subspace are reflected in their proof that this class of permutations is invariant under 
the shift 7Ti7T2 • ■ • 7T n i — > (tti + 1)(7T2 + !)■•■ (7r n + 1), where the addition is modulo n. 



8 Concluding remarks 

It is straightforward to design a Viennot "pyramid" to compute the number a n of S'-avoiding permu- 
tations. For the original Viennot triangle, see [201 ED- Let the entr y a«'"" !m o f th e pyramid be the 
number of permutations in the symmetric group on n elements, avoiding the set S and ending with 
the m entries i±, . . . , i m . Then the entry an''"' 1 " 1 is a sum of entries of the form a > 1 ^y i "' %m ~ 1 . This sum 
is a discrete analogue of the operator T. How far does this analogue between the discrete model and 
the continuous one go? Does the function f n = j 1 ™-"^!) approximate the n-th level of the pyramid? 
More exactly, how well does the integer Un'"' %m compare with n\ ■ f n (ii/n, . . . ,i m /n)7 
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In the case of descent pattern avoidance, can one prove that T restricted to the invariant subspace V 
is compact? We have done so in the case of 123-avoiding permutations. 

Consider the graph G$ of overlapping permutations on the vertex set & m . What is the smallest 
number of edges one has to remove in order to make the graph not strongly connected? Clearly one 
can remove m edges disconnecting the vertex 12 • • • m. Is m the right answer? This would suggest 
that one can remove m — 1 edges without making the directed graph disconnected. 

A more general enumeration problem is as follows. For a function w on 6 m _|_i define the weight of 
a permutation ir = {tt\ti2 ■ ■ ■ vr n ) by the product 

n—m 

Wt(7r) = jj w(II(7rfc,... ,7Tfc +m )). 
k=l 

Now what can be said about the values and asymptotics of the sum 

otiM = ^2 wt ( 7r ) 

TT&6„ 

as n tends to infinity. For instance, if w is a positive function we know by the result of Kreixi and 
Rutman that 

a„(w) ~ c • A n • nl. 

An operator is called a Volterra operator if its spectral radius is 0. What can be said about Volterra 
operators of the form (II. ID ? More specifically, what can be said about the asymptotic of (T n (l), 1)? 
See Examples II. 13[ fl. 141 and 13.91 for different behaviors. 

Are there examples of forbidden sets S such that the associated operator T has non-simple non- 
zero eigenvalues. This situation can also be analyzed using Theorem 22 in [61 Section VII. 3]. When 
the operator T satisfies a symmetry condition as in Lemma 14.71 or in Lemma 14.8} we conjecture that 
all the eigenvalues have index 1. 

In the case of 213-avoiding permutations are all the eigenvalues simple? In other words, are there 
a common non-zero root to the two equations (|6.3[) and (|6.10p ? 

In this paper our object is to understand consecutive patterns avoidance. Generalized pattern 
avoidance was introduced by Babson and Steingrimsson [2j. Is there an analytic approach to obtain 
asymptotics for these classes of permutations? Lastly, it would be daring to ask for an analytic proof 
of the former Stanley- Wilf conjecture, recently proved in [16j . 
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A Spectral Approach to Consecutive Pattern-Avoiding Permutations 



Richard Ehrenborg, Sergey Kitaev and Peter Perry 



Abstract 

We consider the problem of enumerating permutations in the symmetric group on n elements 
which avoid a given set of consecutive patterns S, and in particular computing asymptotics as n 
tends to infinity. We develop a general method which solves this enumeration problem using the 
spectral theory of integral operators on L 2 ([0, l] m ), where the patterns in S have length m + 1. 
Krem and Rutman's generalization of the Perron-Frobenius theory of non-negative matrices plays 
a central role. Our methods give detailed asymptotic expansions and allow for explicit computation 
of leading terms in many cases. As a corollary to our results, we settle a conjecture of Warlimont 
on asymptotics for the number of permutations avoiding a consecutive pattern. 

1 Introduction 

In this paper, we study integral operators of the form 

(Tf)(xi,...,x m )= x(t,xx, . . . ,x m ) • f(t,xi, . . . ,x m -i)dt (1.1) 
Jo 

and their application to enumerating permutations that avoid a consecutive pattern. Here x is a real- 
valued function on [0, l] m which takes values in [0, 1] and is continuous away from a set of measure 
zero in [0, l] m+1 . As we will show, operators of this type arise naturally when counting permutations 
that avoid a consecutive pattern of length m + 1. 

To define the enumeration problem, let 6 n denote the symmetric group on tl elements. For tt £ &n 
we write tt = (ttitt2 ■ ■ ■ ir n ) where the tt^ are the integers from 1 to n. For x = (xi, . . . , x n ) S R n with 
Xi ^ Xj for all distinct i and j, we denote by H(x) the unique permutation tt £ & n with tt; l < 7Tj if and 
only if xi < Xj. A pattern of length k is an element a of If it is a permutation of length n > k, it 
avoids a if n(7Tj, ftj+i, • • • , Kj+k-i) 7^ & for all j with \ < j < n — k + 1. More generally, if S C <& k we 
say that tt avoids S if it avoids each a £ S. That is, S is the set of forbidden patterns. 

For a subset S of <5 m +i, denote by a n (S) the number of permutations tt £ & n that avoid S. 
Observe that for n < m we have a n (S) = n\. Our goal is to compute asymptotics of a n (S) as n tends 
to infinity. Throughout the paper we will assume that m > 2. 

For S C 6 m+ i we define a function \ on [0, l] m+1 by setting x( x ) = if Xj = Xj for some distinct 
indices i and j, and otherwise setting 

= { i n(x) i s, 
x[x> \ o n(x) e 5. 1 j 

Finally, let T be the operator of the form (II. ip . The standard inner product on L 2 ([0, l] m ), is defined 

by 



(/,5) = / f(x) ■ g(x)dx, 
'[o,i] m 
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and hence the I? norm is given by ||/|| = -J (f, /). Moreover, let 1 denote the constant function 1 
on [0, l] m . It is straightforward to prove (see Proposition 12, 10] ) that the formula 

^1 = (T n - m (l),l) (1.3) 

holds for n > m. Note that the left-hand side is the probability of selecting a permutation ir S & n at 
random that avoids the set S. 

The asymptotic behavior of powers of a bounded linear operator is determined by its spectrum. 
Recall that if A is a bounded linear operator from a Hilbert space to itself, the resolvent set of A is the 
set of all z S C so that (zl — A)' 1 is also a bounded operator. The complement in C of the resolvent 
set is the spectrum of A, denoted cr(A), and the spectral radius of A is given by 

r(^) = sup{|A| :Xea(A)}. 

The peripheral spectrum of A is the intersection of cr(A) and the circle of radius r(A) in the complex 
plane. As we will see, the peripheral spectrum of the operator T consists at least of a real eigenvalue 
at r(A), although this need not be the only eigenvalue in the peripheral spectrum. Finally, define the 
adjoint of an operator A to be the operator A* that satisfies (/, A*(g)) = (A(f),g). 
Using spectral theory, we obtain: 

Theorem 1.1. Let S be a set of forbidden patterns in G m+ i. Then the nonzero spectrum of the 
associated operator T consists of discrete eigenvalues of finite multiplicity which may accumulate only 
at 0. Furthermore, let r be a positive real number such that there is no eigenvalue ofT with modulus r 
and let Ai, . . . , A& be the eigenvalues ofT greater in modulus than r. Assume that Ai, . . . , A& are simple 
eigenvalues, with associated eigenf unctions ipt and that the adjoint operator T* has eigenf unctions ipi 
corresponding the eigenvalues Aj. Then we have the expansion 

a n {S)/n\ = (T"" m (l), 1) = V fo' ^-Mi) . X n-m + Q{r ny (L4) 

{<Pi,V>i) 

Moreover, when the operator T has a positive spectral radius, that is, r{T) > then the spectral radius 
is an eigenvalue of the operator T. 

Observe that when T has spectral radius 0, this result is not useful. 

In applications the error term in equation f|l .4[) can be selected to be the modulus of the next 
eigenvalue, that is, 

k 

a n (S)/n\ = • A? + 0(|A fc+1 |"), (1.5) 
i=i 

where the eigenvalues satisfy |Ai| > • • • > |A&| > |A&+i| and A^+i is a simple eigenvalue. This follows 
from letting r be smaller than |Afc+i| given a few more terms in equation (II. 4j) and increasing the error 
term. 

For certain sets of patterns, we can show that there is a unique largest eigenvalue with respect 
to modulus and that this eigenvalue is simple and positive. Thus we obtain the leading term of the 
asymptotics of a n {S)/n\. A sufficient condition for the peripheral spectrum of T to consist only of 
the simple eigenvalue r(T) is as follows. If (X, (j) is a measure space and / E L 2 (X, //) is a real-valued 
function, we will say that / > if f{x) > for almost every x £ X, and / > if f(x) > for almost 
every x. A bounded operator A on L 2 (X,/j.) is positivity improving if for any / > different from 
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there is an integer k (possibly depending on /) so that A k f > 0. Krein and Rutman |16[ Theorem 6.3] 
showed that, for such an operator, the spectral radius r(A) is a simple eigenvalue r(A) and all other 
eigenvalues A are smaller than r(A), that is, |A| < r{A). Furthermore, the associated eigenfunction ip 
is positive almost everywhere. Finally, the adjoint operator A* also has a largest real, positive and 
simple eigenvalue at r(A) and an almost everywhere positive eigenfunction ip. We summarize this 
discussion in: 

Theorem 1.2. If the operator A is positivity improving then its largest eigenvalue is real, positive and 
simple. 

A sufficient condition for T to be positivity improving can be formulated in combinatorial terms as 
follows. If 5 is a set of patterns, let Gs be the graph with vertex set & m and a directed edge from tt to 
a if there is a permutation r £ 6 m+ i\5 with n(ri, . . . , r m ) = tt and II(t2, . . . , r m+ i) = a. The graph 
G$ is known as the graph of overlapping permutations. Recall that a graph G is strongly connected if 
any vertex of G is connected to any other vertex by a directed path. 

Theorem 1.3. Suppose that the graph Gs is strongly connected and that the monotone permutations 
12 • • • (m + 1) and (m + 1) • • • 21 do not belong to S. Then the operator T is positivity improving. 

As a corollary we have the following result, proving a conjecture of Warlimont [23]. See also 
Theorem 4.1 in [9]. 

Corollary 1.4. Let S consist of a single permutation a. Then the asymptotics of a n (S) is given by 
a n (S)/n\ = c ■ X n + 0(r n ), where c, A and r are positive constants such that A > r. 

Proof. When the permutation a differs from the two monotone permutations 12 • • • {m + 1) and (m + 
1) • • • 21, Theorem 11.31 applies to the forbidden set S = {cr}. The two remaining cases are equivalent, 
and can be settled by Theorem 11.71 □ 

Another application of Theorem 11.31 is the following. Call a permutation tt in © n decomposable 
if there exists an index i such that 1 < i < n — 1 and 7r(l),...,7r(i) < i (equivalently i + 1 < 7r(i + 
1), . . . ,7r(n)). A permutation that is not decomposable is called indecomposable. These permutations 
are also known under the terms connected and irreducible. 

Theorem 1.5. Let S be a subset of G m+ i such that each permutation in S is indecomposable and the 
monotone permutation (m + 1) • • • 21 does not belong to S. Then a n (S) has the asymptotic expression 
a n (S)/n\ = c • A™ + 0(r n ), where c, A and r are positive constants such that A > r. 

More generally we can characterize the spectrum of T in terms of another graph associated to S. 
To define it, let A n for tt £ & m denote the set of points x = (x±, . . . , x m ) £ (0, l) m with Xi / Xj for 
i ^ j and II(x) = ir. The graph H$ has vertex set U 7T ^Q m A n and directed edges from {x\, . . . ,x m ) 
to (x2, • • • , x m+ i) if x\ ^ x m+ i and TL(xx, . . . , x rn +\) ^ S. We show that if the graph H$ is strongly 
connected there is an upper bound on the length of the directed path connecting any two vertices. 
We define the period of a strongly connected graph G as follows. Fix a vertex v of G and, for k a 
non-negative integer, let be the set of all vertices in G that can be reached from v in exactly k 
steps. The set Q of all k with v £ is a semigroup and generates a subgroup dZ of Z. The integer d 
is the period of the graph G. Note that, if G is strongly connected, then G has period d for some 
positive integer d. Finally, a graph is called ergodic if it is strongly connected and has period 1. 
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Theorem 1.6. Suppose that S is a set of forbidden patterns and that the graph Hg is strongly con- 
nected with period d. Then the operator T has positive spectral radius r(T) and T has a simple 
eigenvalue A = r(T) with strictly positive eigenf unction. Moreover, the spectrum of T is invariant 
under multiplication by exp(2 • tt ■ i/d). 

In case when the period is 1, we have the conclusion: 

Theorem 1.7. Suppose that S is a set of forbidden patterns such that the graph Hs is ergodic. Then 
the operator T is positivity improving. That is, the operator has a unique largest eigenvalue which is 
simple, real and positive and the associated eigenfunction is positive. 

For certain explicit patterns, we can compute the spectrum and eigenfunctions of T and obtain 
sharp asymptotic formulas for a n (S)/n\. 

Example 1.8. When S is empty, directly we have ct n (S) = n\ for all n > 0. The associated operator 
only has one non-zero eigenvalue, namely 1, with eigenfunction and adjoint eigenfunction 1. This is 
the only case we know where the number of eigenvalues is finite. 

Example 1.9. If S = {123} we show that the operator T has a trivial kernel and spectrum given by 
i x k} k& where 

A ^ 
k 2.*.(* + i)- 

Furthermore all the eigenvalues are simple. We also compute the eigenfunctions of T and the adjoint 
operator T* and obtain 

a„(123) ( 1 \ yll+1 



nl «P^J-*r+°(|A-lH CM) 

where Ao = r(T) and A_i is the next largest eigenvalue in modulus. For more terms in the asymptotic 
expansion see Theorem 15.41 

Example 1.10. If S = {213}, we show that the nonzero eigenvalues of the operator T are the roots 
of the equation 

erf ivh) = V? 

which has the unique real root Ao = 0.7839769312 .... Moreover, Ao is the largest root in modulus of 
the equation. We then have 



a„(213) 



exp 



n! ~' r V 2 • Xl 



Xq +1 + O (|Ai| n ) (1.7) 



where Ai^ = 0.2141426360 . . . ± 0.2085807022 . . . • i are the next two largest roots of the eigenvalue 
equation. See Section [6] for the calculations. 

Example 1.11. If S = {123, 321}, the numbers a n (S) are given by a n (S) = 2E n for n > 2 where E n 
is the nth Euler number. In this case, we can use spectral methods to obtain the classical convergent 
expansion 

§=2. y, (^y 1 - <L8 » 

j odd 
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This formula was derived by Ehrenborg, Levin, and Readdy [8j (Corollary 4.2) by using Fourier series. 
In this case the spectrum of the operator T is real and invariant under the reflection A h-> —A; in par- 
ticular the two largest eigenvalues are ±2/tt. Moreover, the matrix U introduced in Proposition 12. 1 ll 

is similar to the matrix ^ for a cyclic permutation of order two. 

Example 1.12. Let S = {123, 231, 312} . In this case the asymptotic expansion converges and we can 
conclude that the number of such permutations in S n is given by n ■ E n —i, a result due to Kitaev and 
Mansour |14| . See Section [3 

It is easy to find examples of patterns S for which p(T) = 0. 

Example 1.13. Let S = {132,231}. An S-avoiding permutation has no peaks (viewed as the graph 
of a function from {1, . . . , n} to itself) and it is easy to see that a n (S) = 2 n ~ l . However, observe that 
the proportion a n (S)/n\ is subexponential. It is straightforward to verify that the operator T has no 
non-zero eigenvalues. Also note that the graph Gs is not strongly connected. 

Example 1.14. Let S = {123,213,231,321}. The directed graph Gs is strongly connected but the 
monotone permutations 123 and 321 are excluded. In this case a n (S) = 2 for all n > 2. 

We close our introduction with a brief overview on the subject of pattern avoidance in permuta- 
tions (for more details we refer to [31 [13]). The "classical" definition of a pattern is slightly different 
than the one provided above. We say that a permutation ir avoids a pattern a if ir does not contain a 
subsequence which is order- isomorphic to a. The study of such patterns originated in theoretical com- 
puter science by Knuth [15]. However, the first systematic study was done by Simon and Schmidt [T9] , 
who completely classified the avoidance of patterns of length three. Since then several hundred papers 
related to the field have been published. 

One of the most important results in the subject is the proof by Marcus and Tardos [J7] of the 
so-called Stanley-Wilf conjecture related to the asymptotic behavior of the number of permutations 
that avoid a given pattern. It states that for any permutation a there exists a constant c (depending 
on a) such that the number of the permutations of length n that avoid a is less than c n . 

In this paper we also study asymptotic behavior of permutations avoiding patterns, but we con- 
sider consecutive patterns, occurrences of which correspond to (contiguous) factors, rather than sub- 
sequences, anywhere in permutations. Simultaneous avoidance of consecutive patterns of length 3 
is studied in [12] using direct combinatorial arguments. Other approaches to study consecutive pat- 
terns were introduced recently which include considering increasing binary trees [10J by Elizalde and 
Noy, symmetric functions |18] by Mendes and Remmel, homological algebra [5] by Dotsenko and 
Knoroshkin, and graphs of pattern overlaps [1] by Avgustinovich and Kitaev (see Chapter 5 of [13] 
entirely dedicated to known methods and results on consecutive patterns). 

In [10] asymptotics for the following consecutive patterns is given: 123, 132, 1342, 1234 and 1243. 
These results are obtained by representation of permutations as increasing binary trees, then using 
symbolic methods followed by solving certain linear differential equations with polynomial coefficients 
to get corresponding exponential generating functions, and, finally, using the following result (see 
Theorem IV. 7] for a discussion). 

Theorem 1.15. If f(z) is analytic at and R is the modulus of a singularity nearest to the origin in 
the sense that 

R = sup{r > : / is analytic in \z\ < r}, 
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then the coefficient f n = [z n ]f(z) satisfies 

limsupl/,,1 1 /™ = 1/R. 

n— >oo 

This differs from our method, as we obtain detailed asymptotic expansions that allows for explicit 
computation of leading coefficients in many cases. As special cases of our results, we get more detailed 
asymptotics for some of the results of Elizalde and Noy [TU] ■ Also note that the associated generating 
function is related to our operator T by the identity 

MS) ■ ^ = 1 + • • • + z™- 1 + z m ■ ((/ - zT)-\l), 1) . 

n>0 

From this identity it is clear that the radius of convergence of the generating function is determined 
by the spectrum of T. 

The underlying principle that makes our method work is that one can pick a permutation in <5 n at 
random with uniform distribution, by picking a point (x\, . . . , x n ) in the unit cube [0, l] n and applying 
the function II. This method has been used in [8] to obtain quadratic inequalities for the descent set 
statistics and in [7] to enumerate alternating 2 by n arrays. 

2 The operator T 

We now begin our study of the operator T. 



2.1 Connection with pattern avoidance 

Recall that S is a collection of forbidden patterns of length m+ 1, that is, S is a subset of 6 m+ i. The 
function x is defined on the unit cube [0, l] m+1 by 



X{x) 

and the operator T on £ 2 ([0, l] m ) by 



1 U(x) $ S, 
U(x) G S, 



(Tf)(xi,...,X m ) = I x{t,Xl, . . . ,X m ) ■ f{t,Xl, . . . ,X m -l)dt. 



For n > m, define Xn on the n-dimensional cube [0, l] n by 



Xn{xi, ■ ■ ■ ,x n ) = Y[ X{Xj,...,X m+ j). (2.1) 



This allows us to express powers of our operator T. If k < m then 

(T k f) (x) = / Xm+k(h,-- -,tk,xi, . . . ,x m ) ■ f(ti, ...,t k ,xi,.. .,x m ^ k )dh ---dtf, 
v ' ifo.il fc 



while if k > m, 



(2.2) 



T k f)(x)= I Xm+k(h, ■ ■ ■ ,tk,xi, . . . ,x m ) ■ f(ti, . . . ,t m )dti ■ ■ ■ dt k (2.3) 

' i[o,il fc 
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Proposition 2.1. The number of S-avoiding permutations in G n for n > m is given by a n (S) 
n\ ■ (T"- m (l),l). 



Proof. Equation (I2.ip allows us to conclude that 

Xn(x) = 

Hence by integrating over the n-dimensional cube 



1 if IT(x) avoids S, 
if IT(x) does not avoid S. 



Xn{X)dx - 

[0,1]" 



since each simplex in the standard triangulation of [0, l] n corresponds to a permutation tt S & n and 
each such simplex has volume 1/n!. Using this observation and the identities (|2.2p and (|2.3p . we can 
rewrite the integral as (T n_m (l), 1). □ 

Lemma 2.2. The operator T is a bounded operator with norm at most 1. 

Proof. Apply the Cauchy-Schwarz inequality to equation (jl.ip in the variable t to obtain 

-1 2 

X(t,2!l, • • -,x m ) ■ f(t,Xx, . . .,x m -i)dt 



< / \x(t,xi, . . . ,x m )\ 2 dt ■ \f(t,xi, . . . , x m _i)\ 2 dt 
Jo Jo 

1 



< / \f(t,xi, . . .,x m -i)\ dt. 
Jo 

Now integrating over the variables x±, . . . , x m , we obtain that ||T(/)|| 2 < ||/|| 2 proving the bound. □ 
Lemma 2.3. The adjoint operator T* is given by 

(T*f)(x)= / x(xi,---,x m ,u)f(x2,...,x m ,u)du. 



Proof. Since x is real- valued, we have that 



(T(f),g) = I / x(t,xi,- ■ ■ ,x m ) ■ f(t,xx, . . . ,x m -i)dt ■ g(xi, . . . ,x m )dxi- ■ ■ dx r 
'[o,i] m Jo 



= / f(t,xx,...,x m -i)dt- x(t,xi, . . . ,x m ) ■ g(xi, . . . ,x m )dx m dtdxi- ■ ■ dx m -i, 
J[o,i] m Jo 

proving the lemma. □ 

The spectrum of an adjoint operator A* is given by conjugate of the spectrum of A. That is, if 
A is an eigenvalue of A, then A is an eigenvalue of A*. However, for our operator T, the complex 
eigenvalues come in conjugate pairs. This is proved using that x is real- valued. 



Lemma 2.4. We have that T(f) = T{f) and T*{f) = T*(f). Hence if X is an eigenvalue of T with 
eigenfunction ip, then A is also eigenvalue of T with eigenfunction Tp. Similarly, if A is an eigenvalue 
of T* with eigenfunction ip, then A is also eigenvalue of T* with eigenfunction i/j. 
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2.2 Eigenvalues and asymptotic expansion 



Since the operator T is bounded, we can use spectral analysis in order to explore the operator T. We 
refer the reader to Dunford and Schwarz [6, Chapter VII] for a more detailed exposition. 

As we defined in the introduction, the resolvent set s(T) is the set of complex numbers z such the 
operator (zl — T) _1 exists as a bounded operator. The spectrum o~(T) of T is the complement of the 
resolvent set ?(T). 

The index of a complex number A is the smallest non-negative integer v such that the equation 
(XI — T) u+1 f = implies (XI — T) u f = for all functions /. Informally speaking, for operators on 
a finite-dimensional vector spaces, the index of an eigenvalue is the size of the largest Jordan block 
associated with that eigenvalue. A point A in the spectrum is called a pole T of order v if the function 
R(z; T) has a pole at A of order v. Theorem 18 in [6, Section VII. 3] states that the order of a pole A 
is equal to its index. 

Define the operator E(X) by the integral 



where C is a positive oriented closed curve in the complex plane only containing the eigenvalue A from 
the spectrum cr(T). It follows from [61 Section VII. 3] that E(\) is a projection. 

Lemma 2.5. The operator T m is compact. 

Proof. The operator T m has the form 



Since X2m is a bounded function we conclude that T m is a Hilbert-Schmidt operator, and hence a 



Using Theorems 5 and 6 in [6l Section VII. 4] we conclude: 

Theorem 2.6. The spectrum of T is at most denumerable and has no point of accumulation in the 
complex plane except possibly 0. Every non-zero number X in a(T) is a pole ofT and has finite positive 
index. For such a number X the projection E(X) has a non-zero finite dimensional range and it is given 
by {/ 6 L 2 ([0, l] m ) : (XI - T) v f = 0} ; where v is the order of the pole X. 

Recall that an eigenvalue A is simple if the range of E(X) is one-dimensional. That is, the eigenvalue 
equation Xtp = Tip has a unique solution up to a scalar multiple and the generalized eigenvalue equation 
Xf = Tf + ip has no solution. 

Lemma 2.7. Let X be a simple eigenvalue of the operator T with associated eigenfunction (p. Let ip 
be the eigenfunction of the adjoint operator T with eigenvalue X. Then the projection E(X) is given by 



Proof. Since the eigenvalue A is simple, the range of the projection E(X) is one-dimensional and 
spanned by the eigenfunction (p. Since the projection is continuous we may assume that it has the 
form E(X)(f) = (f,a) ■ (p for some function a. It is straightforward to observe that since E(X) is a 
projection, that is, -E(A) 2 = E(X), we have that (<p,ac) = 1. 





compact operator. 



□ 
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Now the adjoint operator E(X)* is given by E(X)*(g) = (g, </?) • a since 

(g, E(X)(f)) = (g, (/, a)-<p) = (a, /) • (g, <p) = ((g, if) ■ a, /) . 

Hence a belongs to the range of the adjoint operator, that is, it is multiple of the eigenfunction ip. In 
fact, we observe that a is given by Up, ip) • tp. □ 

Proof of Theorem By analytic functional calculus we can evaluate the operator T n ~ m by inte- 
grating in the complex plane; see Theorem 6(c) in [61 Section VII. 3]. We have 



1 



2vri T\ z \=r zl — T 



dz, 



where R is greater than the spectral radius of T and we orient the circle in positive orientation. 

Let a be the set {Ai, . . . , A&} and let E(a) denotes the sum of the projections E{X\) + • • • + E(Xk)- 
By Theorem 22 in [6j Section VII. 3] and that the eigenvalues Ai, . . . , are simple, we have that 



T n - m ■ E(a) = J2 E(Xi) ■ A™"™ 



i=l 



We can estimate the operator T n rn ■ (I — E(a)) by shrinking the path of integration to a circle of 
radius r 

T ". (I-E{a)) = — ■ 4) — —(I:. 



2m 



We bound this integral by 

\\T n ~ m ■ {I - E{o))\\ 



1 




2tH ' 


J\z 



\z\=r 



< 

~ 2vr 



\z\=r 



z n- 


-m 


-_ r zl - 


- T 


z n- 


-m 


-. r zl - 




- T 


1 




zl - 


T 



dz 

dz ■ r n - m 



-i 



< sup (zl — T) 

\z\=r 

= 0{r n ), 

where the last equality follows from that the supremum does not depend on n. Hence the inner 
product (T n_m • (I — E(a))l, 1) is also bounded by 0(r n ). Thus we conclude that 

(T"- m l,l) = (T n - m E(a)l,l) + (T n ~ m - [I - E(a))l,l) 

k 



Y J {E{X l )l,l)-Xr m + 0[r n ) 

+ 0(r n ). 



i=l 
k 



i=l 



□ 



9 



If the eigenvalue A is not real then A is also an eigenvalue. The two terms corresponding to these 
two eigenvalues can be combined as follows. Let A = r • e 1 ' 8 and (ip, 1) • (l, ip) / up, tp) = s ■ e 1 '^ . Then 
we have 

A „_ m | (p,l)-(l,ip) yn- m = 2 .ft( a . e *'P . r n-m , e i-{n~m)-e 

(<p,W) (WiVO ^ 

= 2-s-r n - m -cos(f3 + (n-m)-6). 
2.3 Bounds on the norm and spectral radius 

Proposition 2.8. Let T be an operator of the form M.l\) and suppose that T has a positive spectral 
radius r(T). Then T has a positive eigenvalue A = r(T) with non-negative eigenfunction ip. Further- 
more, the eigenfunction ip is almost everywhere positive. Similarly, the adjoint operator T* also has 
A = r(T) as an eigenvalue with an almost everywhere positive eigenfunction ip. 

Proof. Let K be the cone of non-negative functions in L 2 ([0, l] m ). Since T preserves this cone and T 
has nonzero spectral radius, it follows from Theorem 6.1 of [16] and the fact that the cone of positive 
functions is self-dual that T and T* both have A = r(T) as an eigenvalue with at least one strictly 
positive eigenfunction. □ 

For an example of an eigenfunction ip corresponding to the largest eigenvalue, taking the value 
on a set of measure 0, see Proposition 17. 21 where 123, 231, 312-avoiding permutations are discussed. 



Lemma 2.9. For k > m the norm of the operator T k is bounded above by y/ a m +fc (S) / (m + k) ! . 

Proof. Applying the Cauchy-Schwarz inequality to equation (|2.3p in the variables t\ through tk, we 
obtain 

2 

Xm+k (t 1 , • • • , tk , X\ , . . . , X m ) • f(ti, . . . ,t m )dt\ 



[0,1] 



< / \Xm+k(tl,. ■ ■ ,t k ,xi, . . . ,x m )\ dt ■ \f(ti,...,t m )\dt 
J[0,l] k J[0A] k 



2 



Xm+k(h, ■ ■ ■ ,t k ,xi, . . . ,x m )dt ■ / \f(tx,...,t m )\ dt. 

[0,l] fc J[0,l] m 

Now integrating the variables x\,...,x m over [0, l] m we have ||T fc (/)|| 2 < a m+k {S)/(m + k)\ 
proving that ||T fc || < y/a m +k(S)/(m + A;)!. □ 

Proposition 2.10. For a nonempty set of forbidden patterns S we have that spectral radius of the 
operator T is less than 1, that is, r{T) < 1. 

Proof. We have a2 m (S) < (2 • m)\ implying that ||T m || < 1. The inequality ||T m ' n || < ||T m || n implies 
that r(T m ) = lim 

n — >oo \\T mn \\ ' n < 1 strictly. The result follows by taking the ?nth root. Q 

Proposition 2.11. Let T be an operator of the form M.l\) and suppose that p = r(T) > 0. Then the 
operator T admits a decomposition of the form 

T = p-U + W 

where UW = WU = 0, U has finite- dimensional range, U maps the interior of the cone of positive 
functions into itself, the operator W has spectral radius less than p, and the eigenvalues of U are roots 
of unity including 1 . 
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Proof. The orthogonal decomposition follows from Theorem 8.1 of [TE] applied to the operator A = 
p~ l -T. □ 

In particular, the leading behavior of powers T n is determined by the spectral radius r(T) and the 
finite-rank operator U. 

Theorem 2.12. Let S be a set of forbidden patterns. Then spectral radius of the operator T is given 
by 

r(T)= lim 

n — s>oo \ n! 



Proof. Suppose first that r(T) = 0. From the inequality |(T n l,l)| < ||T n || we immediately conclude 
that (a n (S')/n!) 1//n tends to as n goes to infinity. If r(T) > 0, then by Proposition 12.11] we have 

^ifl = r (T) n OJ n l, 1) + (W n l, 1) 
n! 

where the second term obeys the estimate 

\(W n l,l)\<(r(T)-e) n 

for some e > and all sufficiently large n. Moreover, (U n l, 1) is periodic in n and strictly positive 
since U maps the interior of the cone of positive functions into itself. Thus (U n l, 1) is both bounded 
above and below by strictly positive constants. It follows that lim, t — >OQ {U n l,l) 1/n = 1 so that 
lim n >ao (a n (5)/n!) 1//n = r(T) as claimed. □ 

This result extends Theorem 4.1 of [9], where consecutive patterns consisting of a single permuta- 
tion were considered. Moreover, Theorem 12.121 characterizes the limit in terms of a spectral quantity 
which can be computed in many cases of interest by solving the eigenvalue problem for the integral 
operator T. 



3 Associated graphs 
3.1 The directed graph Hs 

In this section we study the spectrum of T using the infinite graph Hs described in the introduction. 
Recall that A n denotes the open subset of (0, l) m with Xi ^ Xj for i 7^ j and Xi < Xj if and only if 
7r(«) < vr(j), and let 

Thus the complement of X consists of those points x with X{ = Xj for at least one pair of distinct 
indices i and j and hence is a set of measure zero. The graph Hs has vertex set X. Recall that the 
directed edges of Hs connected points x and y in X with Xj + \ = yj for 1 < j < m — 1, x\ ^ y m , 
and II(:ei, . . . , x m , y m ) £ S. It follows from the definition that x(t, x±, . . . , x m ) = 1 if and only if there 
is a directed edge from (t, x±, . . . , x m -i) to (xi, . . . , x m ). That is, the function \ encodes the edge 
information of the graph Hs ■ 

The next lemma connects the graph Hs to mapping properties of the operator T. It will be used 
to show that the operator is positivity improving. 
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Lemma 3.1. Suppose that x,y G X and that there is a directed path from x to y of length k > m. 
Suppose further that f is a non-negative continuous function such that f is non-zero in a neighborhood 
ofx. Then (T k f)(y) > 0. 

Proof. Assume that the directed path is 

x = (xi, . . -,x m ) — ► (x 2 , ■ ■ -,x m+ i) — > ■ ■ • — ► (x k+1 , . . .,x k+m ) = y. 

Let e be the minimum of the following finite set 

{\xi — Xj\ : 1 < i < j < k + m, j — i < m} U {xj, 1— Xi : 1 < i < k + m} . 

Observe that e > by the definition of X. Let 5 = e/3. For Sj G [xi — S, x^ + <5], 1 < i < k + m, we 
have that 

(si,...,s m ) — > (s 2) ... 

is also a directed path in Hs- It follows that Xk+m(si, ■ ■ ■ , Sk+ m ) = 1 for all such s. Using (|2.3p we 
may estimate 

"xi+6 rxk+S 



. . rxi+6 rx k +6 

(T k f) (y) > / ••• / f(t l ,...,t m )dt 1 ---dt k 



= (2S) k - m / ••• / f(t 1 ,...,t m )dt 1 ---dt m 

Jxi—5 Jxm — S 

> o, 

where in the last step we have used the positivity of / in & neighborhood, of (3^1, ■ ■ • ; %m 
Proposition 3.2. Suppose that Hs is strongly connected with period d. Then there is a decomposition 

d-l 



X=[jY i 



i=0 

of X into disjoint sets Y{ with the property that T : L 2 (Yi) — > L 2 (Yi + \), where Yd = Yq. 

Proof. Pick a base vertex v of Hs- Let X k be the set of all vertices in Hs that can be reached from v 
in k steps, and let Q be the subset of the non- negative integers defined by 

Q = {k : v£X k }. 

Then Q is a semigroup under addition and generates a subgroup of the integers Z. A subgroup of Z 
has the form dL for some positive integer d; in this case, d is the period of the graph Hs- 
Now define 

Y t = |J Xj. 

j: j=i mod d 

Observe that every directed edge in the graph Hs goes from some Yi to the next Y,- L+ i (with addition 
modulo d). Also, observe that the sets Y{ are pairwise disjoint. 

We claim that each Yi is open. To see this, suppose that y E Y^. Pick a path from some vertex x 
to the vertex y having length greater than 2m. We can perturb this path in a small neighborhood of 1 
y using a variant of the argument used at the beginning of Lemma IBTTl and conclude that Yi is open. 
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Now pick a permutation it € & m - Since the sets Yq,..., l^-i are pairwise disjoint we have that 
A n = A n n X = U^~ 1 (A 7r (~)Yi). Note A^ is a connected set and A n n 3^ are all open. A connected 
set can only be the disjoint union of one open set and hence there exists a unique index i such that 
Ajr C Y{. Hence each set Yi is the disjoint union of the sets A n . 

Finally, suppose that / is a continuous function on [0, l] m with support in the set Y^. We claim 
that Tf is supported in the next set Yi + ±. To see this, note that x(i,xi, . . . ,x m ) = 1 if and only if 
there is a directed edge from (t, x±, . . . , x m _i) to [x\, . . . , x m ). Assuming that f(t, x\, . . . , x m -i) is 
non-zero where (t, xi, • • • , x m —i ) belongs to the set Yi. Since / is continuous function, / is supported 
in a neighborhood of the point (t, xi, . . . , x m _i). By applying the definition of the operator T, we have 
the function Tf supported in a neighborhood of the point (xi, . . . , x m ) in the next set Yi + \. That is, 
Tf is supported in Yi + \. □ 

The next lemma is a straightforward consequence of the definition of the graph H$ and hence the 
proof is omitted. 

Lemma 3.3. Suppose that a : (0, 1) — > (0, 1) is a strictly increasing function. If 

(xi, . . . , X rn ) y (X2, • • • , Xm+l) > • • • y (Xk+l, • • • , Xk+rn) 

is a directed path in Hs, then so is 

(a(xi), . . . , a(x m )) — y (a(x 2 ), a(x m+1 )) — > ■■■ — > (a(x k+1 ), a{x k+m )) . 

We next show that there is an upper bound on the length of the directed path between any two 
vertices in the case when the graph Hs is strongly connected. 

Proposition 3.4. Suppose that Hs is strongly connected with period d. Then there is a positive 
integer N , a multiple of d, such that for any two points x and z in the same component Yi there is a 
path from x to z in the graph Hs of length N . Especially, between any two vertices in the graph Hs 
there is a directed path of length at most N + d— 1 . 

Proof. For two vertices x and y of the graph H, let D(x,y) denote the length of the shortest path 
from x to y. Observe that this is not a distance since it is not symmetrical in general. 

For each permutation ir £ & m pick a point x w in A^-, such that all its coordinates are greater 
than 1/2. Similarly, for all < i < d — 1 pick a point yi in Y^, such that all its coordinates are less 
than 1/2. Let K denote the maximum of the finite set 

{D(x n , yi) : A n CYi} U {D{y h x T ) : A T C Yi}. 

Note that K is a multiple of the period d. Let Qi denote the semigroup 

Qi = {k : there is a path from to of exactly length k} . 

Note that Qi C d ■ N. Furthermore, there is no multiple e of the period d such that Qi C e • N. Hence 
there exists a positive integer rrii such that d ■ N + mi C Qi. That is, there is path from yi to yi of 
any length which is a multiple of d and greater than or equal to mj. Let M be the maximum of mo 
through ma-i- Note again that M is a multiple of d. 

Let N = 2 • K + M. Pick two permutations ir and r such that A^, A r C Yi. We claim that there 
is a path from the point x n to the point x T of length N. We find this path by picking three paths. 
First, choose a path p\ from x^ to yi of length at most K. Second, choose a path p^ from to x T of 
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length at most K. Finally, pick a path p 2 from yi to of length M + 2 • K — D(x 7r , yi) — D(y^,x T ) 
which is at least M. By concatenating the three paths pi, p 2 and P3, the result follows. 

Now let show that there is a path from any point x in Yi to any other point z in Yi of length N. 
Let 7r and r be the permutations ir = H(x) and r = n(z). We do so using that we did choose 
the coordinates of x n and x T greater than the coordinates of y^. Namely, let e > be the smallest 
coordinate of the two points x and z. Let y be the point £ ■ yi- Now there is a monotone function a 
such that a(x 7T ) = x and a(yi) = y. Similarly, there is a monotone function j3 such that = y 

and /3(x T ) = z. Concatenating the three paths a(pi), e -pi and /3(ps) we obtain a path from x to z of 
length iV. □ 

Proof of TheoremUJh It follows from Proposition that T d : L 2 (Yi) — > L 2 (Yi). We will denote 
by A the restriction of T d to L 2 (lo)- Choosing a positive integer p with p ■ d > m we see that 
A p is a compact operator from L 2 (Yq) to itself. The operator A p has discrete spectrum which may 
accumulate only at so the same is true for A. We claim that, if tp is an eigenfunction of A with 
nonzero eigenvalue X d , then T*<£> is a nonzero eigenfunction of T d with eigenvalue X d and support in Y^. 
The only nontrivial part of this claim is that T l ip ^ 0. To see this, note that T d ~ l (T l ip) = X d ■ ip so 
T*(/? cannot be zero since A 7^ and tp is a nontrivial eigenfunction. 

Suppose now that A d G cr(A), let 93 be an eigenfunction of ^4 corresponding to the eigenvalue A , 
let uj be a dth root of unity, and let 



The function ip is nonzero because the right-hand terms are nonzero and have disjoint supports. A 
direct computation shows that Tip = oj ■ A • tp so the spectrum of T contains all of the numbers A • ui 
where X d G o~{A). On the other hand, any eigenvalue /i of T gives rise to an eigenvalue \i d of T d , so 
the nonzero spectrum of T consists exactly of the numbers oj • A where a; is a dth. root of unity and X d 
is an eigenvalue of T d . 

By combining Lemma 13.11 and Proposition 13.41 we have that for a non-negative, but non-zero, 
function / in L 2 (Y ) that T N (f) = A N l d {f) is a positive function. Hence the operator A is positivity 
improving and by Krem and Rutman |16|. Theorem 6.3] the operator A has a positive spectral radius 
r(A). Furthermore, the spectral radius r(A) is a simple eigenvalue, all other eigenvalues are smaller 
in modulus and the eigenfunction p corresponding to r{A) is positive. Hence the spectral radius of T, 
r(T) = ^/r(A), is positive and a simple eigenvalue of T. Finally, the positive function tp is also an 
eigenfunction of T corresponding to the eigenvalue r(T). □ 

3.2 The directed graph Gs 

We examine relations between the infinite graph Hs and the finite graph Gs- We need the following 
lemmas, 

Lemma 3.5. Let x = (x±, . . . ,x m ) be a vertex in Hs such that U(xi, . . . ,x m ) = ir. Furthermore 
assume there is an edge in Gs labeled r leaving the vertex it. Then there exists x m+ \ in (0, 1) such 
that n(xi, . . . ,x m +i) = r. That is, there is an edge in Hs leaving the vertex x. 

Proof. Observe that r(m+l) is bigger than exactly r(m + l) — 1 of the numbers r(l), . . . , r(m). Hence 
pick x m +i such that it is bigger than exactly r(m + 1) — 1 of the numbers x\, — , x m . 

Iterating this lemma we obtain: 
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Lemma 3.6. Given a directed path from ir to a in the graph Gs- Let x be a vertex in Hg such that 
U(x) = tt. Then this directed path can be lifted to a directed path in Hs that ends with a vertex y such 
that II(y) = a. 

Now we give a sufficient condition for Hs to be strongly connected in terms of the graph Gs ■ 

Proposition 3.7. Let S C @ m+ i ; suppose that Gs is strongly connected, and suppose that the two 
monotone permutations 12 • • • (m + 1) and (m + 1) • • • 21 do not belong to the set S. Then the graph Hs 
is ergodic. 

Proof. We first prove that Hs is strongly connected. Let x = [x\, . . . , x m ) and y = (yi, . . . , y m ) be two 
vertices of Hs, and let tt = H(x) and a = H(y). Since Gs is strongly connected we can find a directed 
path from n to m ■ ■ ■ 21 in Gs- This directed path lifts to a directed path from x to z = (z±, . . . , z rn ) 
in where z\> z 2 > • • • > z m . 

For any directed graph G define the reverse graph G* to the graph G where we reverse the 
direction of each edge. In the reverse graph G* s we have a directed path from a to 12 • • • m. Lifting 
this directed path to a directed path in H* s and then reversing the path, we obtain a directed path 
from w = [w\, . . . , w m ) to y in Hs, where W\ < ■ ■ ■ < w m . 

Finally, there is a directed path from m ■ ■ ■ 21 to 12 • • • m in Gs- Hence we know that there is a 
directed path from u = (ui, . . . , u m ) to v = (v±, . . . , v m ) in Hs, where u± > • • • > u m and v\ < • • • < v m . 

Choose a so that < a < mm(z m , wi). We then have a directed path from a ■ u to a ■ v. Now, 
there is a directed path from z to a • u of length m, namely 

z = Oi, ...,z m ) — > (z 2 , ■ ■ .,z m ,aui) — > ■■■ — > (z m ,aui, . . . ,an m _i) — > (aui, . . .,au m ) = a-u 

using the fact that (m + 1) ■ ■ ■ 21 is not forbidden. We can now concatenate these five directed paths 
to obtain a path from x via z, via a • u, via a • v, via w, to y. 

Since Gs is strongly connected and has ml vertices, an upper bound on the length of the directed 
paths in Gs chosen above is ml — 1. Hence, the path that we have constructed has length at most 
3(m! - l) + 2m. 

To observe that Hs has period 1 note that we can construct a directed path from the vertex x 
to the vertex y that has length one more than the above construction. Namely, the path from tt to 
m ■ ■ ■ 21 can be extended by adding the loop (m + 1) • • • 21 at the end. Now by concatenating these 
two paths with a path from y to x we obtain two cycles whose lengths differ by one. Since the greatest 
common divisor of two consecutive integers is one, the graph Hs is ergodic. □ 

Recall that a permutation tt in 6 n is indecomposable if there is no index i such that 1 < i < n — 1 
and 7r(l), . . . ,ir(i) < i. Otherwise the permutation is decomposable. 

Proposition 3.8. Let S C (3 m+ i such that each permutation r G S is indecomposable. Then the 
graph Gs is strongly connected. 

Proof. Given two vertices tt = (tt\, . . . , TT m ) and a = (a±, . . . , a m ) of Gs- Since the integers tt\, . . . , TT m , 
a\ + m, . . . , a m + m are distinct, the following path (described by its edges) is well defined: 

n(7ri, . . . ,TT m ,o-\ + m),n(7r 2 , . . . ,TT m ,ai + m,a 2 + m), . . . ,n(7T m ,cJi + m, . . . ,a m + m). 

This path goes from the vertex tt to vertex a. Note that every edge r' on the path is decomposable. 
Hence this path avoids the forbidden indecomposable edges of S. □ 
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Proof of Theorem \1.5[ This follows directly from Theorems 11.21 and 11.71 and Propositions 13.71 and 13.81 

□ 

Proof of Corollary \l-4\ Let r be the single permutation in the set S. If t is one of the two monotone 
permutations the result will follow from descent pattern avoidance, see Theorem 14.21 If r is indecom- 
posable the result follows from Theorem 11.51 Finally if r is decomposable apply Theorem 11.51 to the 
upside down permutation {m + 2 — r(l), . . . , m + 2 — r(m + 1)) which is not decomposable. □ 

Note that there are examples of patterns S so that H$ does not have cycles even though the 
graph Gs has cycles. 

Example 3.9. S = {312,321}. In this case the graph Gs has a cycle. However, the graph H$ 
does not have a cycle and hence is not strongly connected. We observe this by noting that for a 
directed edge (x,y) — > (y,z) in Hs we have that x < max(y,z). Hence none of Xj's in a /c-cycle 
(x\, X2) — > (x2,X3) — > ■ ■ ■ — > (xk, x\) — > (xi,X2) can be the largest. 

Via the classical bijection tt 1 — > tt (see [20^ Section 1.3]) one obtains that the number {312,321}- 
avoiding permutations are in bijection with involutions, that is, permutations ir such that tt 2 = id. 
This was first observed by Claesson [3J. It follows that the generating function is exp(z + z 2 /2) and 
the asymptotic is 1/V2- exp(— 1/4) • (n/e) n / 2 • exp(y / n). 

On the other hand, if Gs does not have a cycle, then neither does H$- 

Lemma 3.10. Let S C © m +i and suppose that Gs has a directed cycle that contains the two vertices 
12 • • • m and m---21. Moreover, assume that the two monotone permutations 12- • • (m + 1) and 
{m + 1) • • • 21 do not belong to the set S. Then the graph Hs contains a directed cycle. 

Proof. Pick a vertex x = (x±, . . . , x m ) such that H(x) = 12 • • • m. The directed path from the vertex 
12 • • • m to the vertex m ■ ■ ■ 21 can be lifted to a path in Hs from the vertex x to a vertex y = 
(yi, . . . , y m ) where Tl(y) = m ■ ■ ■ 21. Similarly, pick a vertex z = (z±, . . . , z m ) such that H(z) = m ■ ■ ■ 21. 
The directed path from the vertex m ■ ■ ■ 21 to the vertex 12 • • • m can be lifted to a path in Hs from 
the vertex z to a vertex w = (w\, . . . , w m ) where IL(w) = 12- • • m. 

Using the two monotone functions a,/3 : (0, 1) — > (0, 1) defined by a(x) = (x + l)/2 and [3(x) = 
x/2, we have two directed paths: one from a(x) to a(y) and one from (3(z) to /3(w). 

Using that (m + 1) • • • 21 is an edge in Gs we have the following path from a(y) to 0(z), namely 

a (y) = {a(yi),---,a(y m )) — > (a(y 2 ), • . . , a(y m ), f3(zi)) — > > (P{zi), ■ ■ ■ , P{z m )) =/?(»• 

Similarly, we have the directed path 

P(w) = (f3(wi), . . . ,P(w m )) — > (f3(w 2 ), ■ ■ ■ ,(3(w m ),a(xi)) — > ■■■ — > (a(xi), . . . ,a(x m )) = a(x). 
Concatenate these four directed paths to obtain a directed cycle in Hs- □ 

4 Computational techniques 

In this section we discuss descent pattern avoidance which is a special case of pattern avoidance. First 
we introduce an analogue of the de Bruijn graph Djj, which has the advantage that it is smaller than 
the graph Gs- Moreover, if the graph Djj is ergodic so is the graph H S (jj) and we obtain that the 
associated operator is positivity improving. Second, for descent pattern avoidance we obtain that the 
eigenfunctions has a simplified form. Finally, in the last subsection we consider pattern avoidance 
that has symmetry. In these cases we show that we can obtain the adjoint eigenfunctions from the 
eigenfunctions. 
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4.1 Descent pattern avoidance 

The descent set of a permutation tt in the symmetric group on n elements is the subset of {1, . . . , n— 1}, 
given by {i : 7Tj > 7Tj+i}. An equivalent notion is the descent word, defined as follows. The descent 
word of the permutation tt is the word u(ir) = u\ ■ ■ -u n -i where U{ = a if tti < 7Tj + i and Uj = b 
otherwise. 

Let U be a collection of aft-words of length m. The permutation tt avoids the set U if there is no 
consecutive subword of the descent word of tt contained in the collection U. 

Descent pattern avoidance is a special case of consecutive pattern avoidance. For instance, per- 
mutations avoiding the word aab is the permutations avoiding the set S = {1243, 1342,2341}, since 
these three permutations are the permutations with descent word aab. More formally, for U a subset 
of {a, b} m define S(U) C & m+1 by 

S(U) = {tt £ & m+1 : u(tt)£U}. 

Then the set of permutations avoiding the descent words in U is the set of permutations avoiding 
S(U). 

For U a subset of {a, b} m define the associated de Bruijn graph Djj by letting the vertex set be 
{a,b} m ~ l . For x,y £ {a, 6} and u E {a, b} m ~ 2 such that xuy U let there be a directed edge from 
xu to uy. When the set U is empty, the graph D\j is the classical de Bruijn graph D m ^\. 

Lemma 4.1. Let U be a subset of {a, b} m . If there is a cycle c of length N > n + 1 that do not consists 
only of the loop a m or not only of the loop b m , then the cycle c can be lifted to a cycle of length N in 
the graph Hgnjy 

Proof. Let c be the cycle 

VlV2 ■ ■ ■ V m -1 ► -"2^3 ■■■V m >-" > V N Vi ■ ■ ■ V m - 2 ► ^2 • • • 

where each Vi is either a or b. We would like to pick N real numbers x±, . . . ,xn in the interval (0, 1) 
such that 

Xi < x i+ i if Vi = a, , . 

Xi > x i+ i if vi = b, 

where all the indices are modulo N. Since all the letters v\ through vn are not the same, we may 
assume that vn-i = a and vn = b. Pick x\ arbitrarily. Pick x 2 through xn-i such that inequality (|4.ip 
is satisfied. Finally, pick x^ in the interval (max(xjv-i,^i)) 1)- Now in the graph Hgnj\ we have the 
cycle 

(xi,X 2 , Xm) > (X2,X 3 , . . . , Xm+l) > ' ' ' > (XN, %l, ■ ■ ■ , X m -l) > (xi,X 2 , ■■■ , X m ). 

□ 

Theorem 4.2. Let U be a subset of {a,b} m . If the de Bruijn graph Djj is strongly connected, then 
the graph Hs(U) * s a ^ so strongly connected. Furthermore, the de Bruijn graph Djj has the same period 
as the graph Hgnn . 

Proof. Note that the graph Djj has the directed edge a m ~ 1 b since otherwise it would not be strongly 
connected. Given two vertices x and y in Hgmy To prove that Hgnj\ is strongly connected it is 
enough to find a directed path from x to y. 
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We can find a path from u(U(x)) to a m_1 in the graph Djj that consists of at least m + 1 edges. 
Similarly to the lifting lemma, Lemma |3.6| we can lift this path to a path in Hs(jj) that starts at the 
vertex x and ends, say, in the vertex z = (z±, . . . , z m ). Note that Z\ < ■ ■ ■ < z m . Moreover, we can 
find a path from a m ~ 2 b to u(H(y)) in Djj that has length at least m + 1. Lift this path to a path that 
ends in the vertex y and begins at w = (w\, . . . , w m ), where w\ < ■ ■ ■ < w m ^\ > w m . 

Let Vi = max(zj, for 2 < i < m. Observe that we have the string of inequalities z\ < V2 < 

■ ■ ■ < v m > w m . We can now concatenate these two paths as follows. Replace each occurrence of Z{ 
and Wi-i by Vi for 2 < i < m in each of the two paths. Then we may connect the vertex (zi,V2, ■ ■ ■ , v m ) 
with the vertex (v2, ■ ■ ■ ,v m , w m ) via the edge that goes across the edge with descent word a m ~ 1 b. Thus 
the graph H S (u) is strongly connected. 

Since there is a graph homomorphism from Hg^jj-j to D\j we know that the period of D\j divides 
the period of H S (jj)- To see that the periods are equal, pick a vertex w of Djj that differs from a m ~ x 
and b m ~ 1 . Then any cycle of length greater than m + 1 through the vertex w in Djj lifts to a cycle 
of the same length in Hg^jy Hence the greatest common divisor of lengths of cycles through w is 
a multiple of the greatest common divisor of the cycle lengths of Hgmy Hence the two periods are 
equal. 

Note that this argument only works when m > 3 since there is no such vertex w in the m = 2 
case. But the remaining m = 2 case is straightforward to check. □ 

4.2 Invariant subspace for descent pattern avoidance 

For an a6-word u of length m — 1 define the descent polytope P u to be the subset of the unit cube [0, l] m 
corresponding to all vectors with descent word u. That is, 

P u = {(xi,. . . ,x m ) £ [0, l] m : Xi < x i+ i if Ui = a and x { > x i+ \ if m = b}. 

Observe that the m-dimensional unit cube is the union of the 2 m_1 descent polytopes P u . Now the 
operator T corresponding to the descent pattern avoidance of the set U has the following form. For 
an a6-word u of length m — 2 and y £ {a, b} we have 

T(f)\ Puv = H X(auy)- f(t,x 1 ,...,x m „ 1 )\ PaK dt (4.2) 



o 



+ / X{buy) ■ f(t,xi,...,x m -i)\ P dt, 



where by abuse of notation we let x( w ) = 1 if w does not belong to the set U and x{ w ) = otherwise. 

Proposition 4.3. Let T be the operator associated with a descent pattern avoidance and k is an 
integer such that < k < m — 1. Let u be an ab-word of length m — 1. Then the function T k (f) 
restricted to the descent polytope P u only depends on the variables x\ through x m _/ t . 

Proof. Proof by induction on k. When k = there is nothing to prove. When 1 < k < m — 1, we know 
by induction that the restriction of T k ~ 1 (f) only depends on x±, . . . , x m _k+i- By the shift of variables 
in the right hand side of equation (|4.2p . we obtain that T k (f) does not depend on the variable x m -k+i, 
completing the induction. □ 

Corollary 4.4. Let T be the operator associated with a descent pattern avoidance and let ip be an 
eigenfunction associated with a non-zero eigenvalue A. Then the eigenfunction restricted to each 
descent polytope P u only depends on the variable x\. 
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Proof. Since X m ■ tp = T m (tp) the eigenfunction has the required form. □ 

Corollary 4.5. Let T be the operator associated with a descent pattern avoidance and let tp be an 
eigenfunction associated with a non-zero eigenvalue A. Assume that f is a generalized eigenfunction, 
that is, it satisfies the equation A • / = T(f) + tp. Then the function f restricted to each descent 
polytope P u only depends on the variable x\ . 

Proof. By induction on k. Assume that 1 < k < m and that / restricted to each descent polytope only 
depends on the variables x\ through x^. Then T(f) + tp only depends on x\ through x^-i showing 
that A • / only depends on x\, . . . , xt-i- □ 

Let V be the subspace of L 2 ([0, l] m ) consisting of all functions / that only depend on the variable x\ 
when restricted to each of the descent polytopes P u . Observe that the subspace V is invariant under 
the operator T. That is, the operator T restricts to the subspace V. Moreover the constant function 1 
belongs to V. Hence to understand the behavior of T n (l) it is enough to study this restricted operator. 

In order to describe the subspace V more explicitly define for an aft-word u of length m—1 the 
polynomial h{u]X\) as follows: 



h(u;xi)= / ldx2 ■ ■ ■ dx m . 

J (xi,X2,..;X m )eP u 

These polynomials were first introduced and studied in [8], with different notation. 

Let p be a vector (p u (%i)) ue { a M.m-i- That is, the vector p consists of one-variable functions in the 
variable x\ and is indexed by ah- words of length m—1. Consider the function / on [0, l] m defined by 

f{x 1 ,...,X m )\p u =p u (x 1 ) 

for all afc-words u of length m—1. Observe that the function / belongs to L 2 ([0, l] m ), and hence to 
the invariant subspace V, if and only if 

l 

h(u;xi) ■ \p u (xi)\ 2 dx\ < oo 



o 



for all ab- words u of length m — 1. For two functions / and g in the subspace V, corresponding to the 
two vectors (p u (xi)) ue ^ b y m -i and (q u (xi)) ue ^ b y m -i, the inner product is given by 



(fi9)= "52 / h(u;xi) -p u (xi) ■ q u (x\)dxi. 



We end this section by a structural result about the subspace V. 

Proposition 4.6. The invariant subspace V is isometrically isomorphic to the Hilbert space 

Proof. The isomorphism of the Hilbert spaces V — > L ([0, 1]) is given by 

(Pu(a?l))« e {o,6}»»-i >• (VHuiX!) -pu(xi)) . 

□ 
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4.3 Symmetries 

Let J and R be the following two involutions on the space L 2 ([0, l] m ): 

(J/)(xi,x 2 , • • • ,x m ) = f(l - x m , . . . , 1 - x 2 , 1 - Xl), 

(Rf){x\,X2, . . . ,X m ) = f(x m ,...,X 2 ,Xi). 

Observe that both J and R are self adjoint operators. 

Lemma 4.7. Assume that x has the symmetry 

x(x\,x 2 , ■ ■ ■ ,x m ,x m+ i) = x(l - x m+ i, 1 - x m , . . . , 1 - x 2 , 1 - xi). 

Then the adjoint of the associated operator T is given byT* = JTJ. Moreover, ifip is an eigenfunction 
of the operator T with eigenvalue A then ip = Jip is an eigenfunction of the adjoint T* with the 
eigenvalue X. Furthermore, we have the equality (l,V0 = (v 9 ;!)- 

Proof. We have that 

JTJf(xi,x 2 ,...,x m ) = JT/(1 - x m , ■ ■ ■ , 1 - x 2 , 1 - xi) 



J / x(t, xi,..., x m ) ■ f(l - x m -i, . . . , 1 - xi, 1 - t)dt 
Jo 

/ x(t, 1 - x m , . . . ,1 - xi) ■ f(x 2 , x m , 1 - t)dt 
Jo 

/ X(l - t, 1 - x m , . . . , 1 - xi) ■ f(x 2 , ...,x m , t)dt 
Jo 



X(xi, ...,x m ,t)- f(x 2 , x m ,t)dt 
T ■ ■ • j x m — i, x m ). 







For the second statement consider the following line of equalities T* Jip = JTJ Jip = JTip = A • Jip. 
Lastly, = (1,7$) = (l,Jlp) = {Jl,Tp) = {l,Tp) = {ip,l). □ 

Similarly to Lemma 14.71 we have the next lemma. Its proof is similar to the previous proof and 
hence omitted. 

Lemma 4.8. Assume that x has the symmetry 

X(XI,X 2 , . . .,X m ,X m+ i) = x(x m +l,X m , ■ ■ ■ ,X 2 ,Xl). 

Then we have that the adjoint of the associated operator T is given by T* = RTR. Moreover, ifip is an 
eigenfunction of the operator T with eigenvalue A then tp = Rip is an eigenfunction of the adjoint T* 
with the eigenvalue A. Furthermore, we have the equality = (v 9 )!)- 

5 123- Avoiding permutations 

A 123-avoiding permutation is a permutation ir G & n with no index j so that ttj < vrj + i < TTj+ 2 , where 
1 < J ' < n — 2. Let a n (123) denote the number of 123-avoiding permutations in <5 n . 
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5.1 Eigenvalues and eigenf unctions 

Since 123-avoiding permutations can be viewed as permutations with no double descents Corollary 14.41 
allows us to recast the problem of finding eigenfunctions in two variables into finding two one-variable 
functions. 

Proposition 5.1. The eigenvalues of the operator T are given by 



V3_ 

2-7T-(fc + i) 



X k = - _ — , lv (5.1) 



where k G Z and the associated eigenfunctions are given by 

x cos(f + f -f) <fO<x<y<l, 

tp k = exp -— - ■ { ) ( (5.2) 



2 • A 7 I sin ( f + ^ • § ) £/ < y < x < 1. 



Proof. Avoiding the pattern 123 is equivalent to avoiding the descent set pattern aa. Hence Corol- 
lary H3J states that the eigenfunctions <p can be written as 



p(x) if < x < y < 1, 
g(jc) ifO<y<x<l. 



Then the defining equations for eigenvalues and eigenfunctions reduces to the integral system: 

X-p{x) = [ q(t)dt, (5.3) 



/ p(t)dt+ [ q(t)dt. (5.4) 
Jo Jx 



A • q(x 

First, differentiating with respect to x, we obtain the first-order system 

X-p(x) = —q(x), (5-5) 
X-q'(x) = p(x) — q(x). (5-6) 

These equations have only the trivial solution if A = 0, so A = is not an eigenvalue. If A ^ then 
the first-order system (I5.5p ~ (l5.6p implies the second-order equation 

A 2 • p"(x) + A • p'{x) + p(x) = 0. 

This equation has the general solution 

oj \ ( 



p(x) = A ■ exp ( — • xj + B ■ exp ( — • x ) , (5.7) 



where u = exp (^-f-^)- That is, ui satisfies the relation oo 2 + to + 1 = 0. Moreover, equation ()5.5 
implies that 



-to ■ A - exp ^— • x*J — to 2 • B ■ exp ( — ■ x ) . (5. 
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Setting x = and x = 1 in equations (|5.3|) and (|5.4|) and using that A/Owe obtain the boundary 
conditions: 

p(0) = g (0), (5.9) 
p(l) = 0. (5.10) 

Substituting the expressions for p(x) and q{x) from equations (|5.5p and (|5.6p into boundary condi- 
tion (|5.9|) we obtain A + S = — cj • A — u 2 ■ B. This is equivalent to u ■ A + B = 0. Hence we may 
set A = 1/2 • exp (^) and £? = A = 1/2 • exp (— Substituting equation ()5.5|) into the second 
boundary condition (|5.10p implies that 



A ■ exp ^— ^ = —5 • exp = u • A - exp 



Cancelling A on both sides and taking the logarithm gives 

UJ 2 • 7T • i UJ 1 

where is an integer. Since cj — w 2 = v3 • i we obtain expression (|5.ip . Moreover p(x) is given by 

1 (n;-i\ ( x\ 1 / 7T-A / 2 x 

= --exp^— j-exp^- J + --exp^-— j-exp^ ■- 

/ X \ / 7T \/3 X \ 

Now equation (|5.5p implies the claimed expression for g(x). □ 

Note that the eigenvalues are ordered by 

A > -A_i > Ai > -A_2 > A 2 > -A_3 > A 3 > • • • > 0. 

Furthermore, the calculations in Proposition 15. II showed that they all have a unique eigenfunction. It 
remains to show that they are simple. 

Proposition 5.2. The eigenvalues A of the operator T are simple, that is, they have index 1. In other 
words, there is no function f(x,y) such that \k ■ f = T(f) + (pi-. 



r(x) if < x < y < 1, 
s(x) if < y < x < 1. 



Sketch of proof. Using Corollary 14.51 we can write the function / as 

f(x,y) = 

Then equation reduces to the integral system: 

A • r{x) = [ s(t)dt + p(x), 
A • s(x) = 



r(t)dt+ I s(t)dt + q(x). 

J X 
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Note that we obtain the two boundary conditions r(l) = and r(0) = s(0). Next differentiating with 
respect to x, we have the first-order system: 

A • r (x) = —s(x)+p'(x), 
X-s'(x) = r(x) — s(x) + q'(x). 

This system of differential equations can be solved as in the proof of Proposition 15.11 However, the 
solution does not satisfy the boundary conditions, completing the sketch. 
By applying the involution J we obtain the adjoint eigenfunction 

cosf^ + ^-i^O if 0<x <y<l, 



2-A / sinf^ + ^-i^) ifO<y<x<l 



Proposition 5.3. Let A be an eigenvalue of T with eigenf unctions (p and let ip be the eigenfunction 
of the adjoint operator T* with eigenvalue A. Then the following identities hold: 

(¥>,!) = (1,?) = ^-A 2 , (5.12) 

faty = |.(-l)*.A.«p(-iV (5.13) 

In particular 



(*>,?) 1 ' ' V2-A 

Proof. In the following calculations we use the facts that cos(\/3/(2 • A)) = (— l) fc /2 and sin(\/3/(2 • 
A)) = (— l) fc \/3/2. We also use the expression for ip in equation (|5.2p . First, we note that 

(cp, 1) = / p(x)dxdy + / q(x)dxdy 

J0<x<y<l J0<y<x<l 
■1 H 

(1 — x)p(x) ■ dx + x • q(x)dx. 
Jo Jo 

Explicit computation shows that 

(l-x)-p(x)dx = ^.A 2 -(l-(-l) fc .exp(-l/(2-A))) 

x-q{x)dx = ^-A 2 -(-l) fc -exp(-l/(2-A)) 
which shows (|5.12[) . Next, using (|5.2p and (|5.1ip and dropping subscripts as before, we have 



(<p,ip) = / p{x) ■ p(l - y)dxdy + q(x) ■ q(l - y)dxdy 

J0<x<y<\ J0<y<x<l 

= J (p(x) - J p(i - y)dy + q{x) - J <?(i - y)dy^j dx. 

Carrying out the y integration and simplifying, we obtain 

{<p,V) = ^-(-l) k - j\-exp (-^rj) dx 
which gives f)5. 13|) . □ 
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5.2 Asymptotics 

The above computations show that all eigenvalues of T are simple and give the eigenvalues and the 
coefficients explicitly. We thus obtain the following expansion for (T n_2 (l),l) = a n (123)/n! as an 
immediate consequence of Theorem 11.14 Propositions 15,11 and 15.31 

Theorem 5.4. Let K be a non-negative integer. The number of 123-avoiding permutations satisfies 
the following asymptotic expansion 



corresponding to the 2K + 1 largest eigenvalues. 

6 2 13- Avoiding permutations 

A 213-avoiding permutation is a permutation tt £ & n which contains no sequence of the form 

TTj +1 < TTj < TTj +2 

for any j with 1 < j < n — 2. We denote the number of 213-avoiding permutations of <3 n by a n (213). 
Thus, S consists of the single permutation 213 and 



By symmetry, the study of 213-avoiding permutations is equivalent to 132-avoiding permutations, 
231-avoiding permutations and 312-avoiding permutations. However the case of 213-avoiding permu- 
tations gives the most straightforward equations. 

6.1 Eigenvalues and eigenfunctions 

In what follows, we will make use of the error function 




where is given 




sum contains 2K + 1 terms 





(6.1) 



which extends to an entire function on C, and the function 




(6.2) 



Let 





Now we characterize the nonzero eigenvalues and eigenfunctions. 
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Proposition 6.1. The non-zero eigenvalues A of the operator T satisfy the equation 

and the corresponding eigenf unctions are 

' Q( x )~J- nq(t) dt ifO<x<y<l, 
q(x) if < y < x < 1, 



where q(x) is given by A6.2\) . 

Proof. The defining relations for the eigenfunctions are 



X-p{x,y) = [ p(t,x)dt+ [ q(t,x)dt, (6.4) 

JO Jy 
rx rl 

X-q(x,y) = / p(t,x)dt+ / q(t,x)dt. (6.5) 

JO Jx 

Now observe that in the right-hand side of equation (|6,5p there is no dependency on the variable y. 
Hence we may replace q(x,y) with q(x). Now subtract equation (]6.5p from equation (|6.4|) 

rv 

A • (p(x, y) - q(x)) = - I q(t)dt. 



That is, 

1 [ y 

p[x,y) = q(x) - - J q(t)dt. (6.6) 



Substitute equation (|6.6p into equation (|6.5D : 



A • q(x) = J (q(t) - i J q{s)ds^j dt + j q(t)dt 



q(t)dt — - • / / q(s)dsdt 
o A J J t 

P 1 PX f k S 

/ q(t)dt - - ■ / / q(s)dtds 
Jo * J J 



q(t)dt — — ■ j s ■ q(s)ds. 



Hence we have the following integral equation for q{x) 

rl rx 

A 2 • q(x) = A • / q(t)dt - s- q(s)ds. (6.7) 
Jo Jo 

Differentiating once we have 

A 2 • q'(x) = —x ■ q(x). (6.8) 
The solution to this differential equation is 

q( x ) = C-e X p(~^\. (6.9) 
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By setting the constant C to be 1 we obtain a solution to equation (|6.8|) . Now substitute this solution 
for q(x) into the integral equation (|6,7p and set x = 0: 



A 2 = A • / exp ( I ill 



o 



2 • A 2 

1/(V2-A) 



v2 • A • / exp (— u 2 ) 

^•A 2 



V2 



erf (1/(V2 • A)) 



where the substitution in the integral is u = t/{y2- A). Hence the non-zero eigenvalues A satisfy 
equation (I6.3p . □ 



Proposition 6.2. All the non-zero eigenvalues of the operator T are simple, that is, there are no 
non-zero eigenvalue whose index is greater than or equal to 2. 

Proof. Assume that A is a non-zero eigenvalue of index at least 2, that is, there is a function f(x,y) 
such that A • / = T(f) + <p. Letting f(x,y) = r(x,y) for < x < y < 1 and f(x,y) = s(x,y) for 
0<y<x<lwe have 

rx pi 

X-r(x,y) = / r(t,x)dt+ / s(t, x)dt + p(x,y), (6.10) 
JO Jy 
px pi 

X-s(x,y) = / r(t,x)dt+ s(t,x)dt + q(x). (6.11) 

JO Jx 

Note that in the right-hand side of equation ()6.11j) there is no dependency on the variable y. Hence 
s(x,y) is a function of x only, and we write s(x) henceforth. Subtracting equation (I6.1ip from equa- 
tion (|6.10p and dividing by A gives 

I rv i 
r(x,y) = S ( X )~J- s{t)dt + -'(p(x, y)-q(x)) 



*(*)- p ■ f (X-s(t)+q(t))dt, 



where the last step is by equation (|6.6p . Substituting this expression into equation (|6.1ip yields the 
integral equation for s(x): 

A • s(x) = J ^s(t) - ^ • (A • s(u) + q(u))du^ dt + J s(t)dt + q(x) 



1 

A 2 
1 

A 2 
1 

A 2 
1 

A ' 



(A • s{u) + q(u))dudt + f s(t)dt + q(x) 
o Jt Jo 

x pu pi 

/ (A • s(u) + q{u))dtdu + I s{t)dt + q(x) 
1 1 Jo Jo 

1 



u ■ (A • s(u) + q{u))du + / s(t)dt + q(x) 

Jo 

[ u-s(u)du+[ s(t)dt-~-[ q(t)dt + 2- q(x), (6.12) 
Jo Jo * Jo 
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where the last step is using equation (|6,7|) . Differentiate with respect to x to obtain the first order 
differential equation in s(x): 

A • s (x) = — — • x ■ s(x) + 2 • q'(x). 
A 

The general solution to this differential equation is 

X I X 



s( x ) = -ja ■ ex P y-JTJi J + C ' 

where we used that q(x) = exp ^— Observe that the homogeneous part of this solution is 

expected. It corresponds to the homogeneous part of the equation A • / = T(f) + <p. Hence we may 
set C = without loss of generality. Setting x = in (|6.12p yields 

0= f s(t)dt-~- [ q(t)dt + 2- q(0). 
Jo A J 

Using that q(0) = 1, q(t)dt = A and 

s(t)dt = — ■ exp ( — - 1 ) - a h- ■ erf 



A 2-AV V2 Vv/2-A 



we obtain the equation 

1 



= — • exp 



A 2 • A 2 

which has no solutions. Hence there is no non-zero eigenvalue of index greater than 2, that is, all 
non-zero eigenvalues are simple. □ 

For completeness we state: 
Lemma 6.3. A function in the kernel of the operator T has the form 



P(x,y) if0<x<y<l, 
ifO < y < x < 1, 



where p(x,y) satisfies Jq p(t,x)dt = 0. 
The adjoint operator T* is given by 

Jo <?(y> u ) du + ly PiVi u ) du if < x < y < 1, 

Jo Q(Vi u ) du + ly u ) du if0<y<x<l. 

Proposition 6.4. For a non-zero eigenvalues A of the operator T the corresponding eigenfunction of 
the adjoint operator T* is 



T*(f(x,y)) 



where 



P*(y) ij0<x<y<l, 
P*(y) ~ J ■ !lv*{u)du ijO<y<x<l, 



p*(y) = -2 • y ■ exp (j^) + 2 • A + • y ■ exp U^) ' erf (^) ■ (6.13) 
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Proof. The defining relations for the eigenfunctions are 

ry 



X-p*(x,y) 
A • q*{x,y) 







*(y,u)du+ / p*(y,u)du, 



*(y,u)du+ / p*(y,u)du. 



(6.14) 
(6.15) 



Observe that there is no dependency on the variable x in equation (|6.14p . Thus we write p*{x,y) 
p*(y). Subtracting these two equations we have 



A- (q* (x,y) -p*(y)) = - / p*{u)du, 

J X 



such that 



f{x,y)=p*(y)---l p*{u)du. 



(6.16) 



Substituting this expression into equation (|6.14p one obtains 



X-p*{y) 



p*{u) — — ■ I p*(v)dv J du+ I p*(u)du 



p*{u)du — — 
A 



1 f y f 1 



Jy 



p* (v)dvdu 



p*(u)du— — -y- I p*(v)dv. 



That is, p*(y) satisfies the integral equation 

X 2 -p*(y) = X- / p*{u)du — y- I p*(v)dv. 
Jo Jy 

Differentiating this equation twice we obtain 

A 2 -y>*'(y) = y-p*(y)- p*(v)dv, 

Jy 

A 2 -P*"(y) = yp*'(y) + 2-p*(y). 
The solution of this differential equation is given by 

,2 



P* (y) = Ci ■ y ■ exp 



y 



2- A 2 



+ c 2 - 



2 • A + V2 • vr • y ■ 



cxp 



2 • A 2 



• erf 



a/2 • A 



(6.17) 

(6.18) 
(6.19) 

(6.20) 



Setting y = in equations (|6.17p and (|6.18p we obtain A • p*(0) = Jq p*(u)du = — A 2 • p'(0). Inserting 
this condition into the solution of the differential equation (16.2Q[) we obtain C\ = — 2 ■ C%. Moreover 
setting C2 = 1 we obtain equation (|6.13p . □ 



Lemma 6.5. A function in the kernel of the adjoint operator T* has the form 

ip(x,y) = 

where q*(x,y) satisfies Jq q*(y, u)du = 0. 



if0<x<y<l, 
Q*(x,y) if0<y<x<l 
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Proposition 6.6. For a non-zero eigenvalue A with eigenfunction ip and eigenfunction tjj of the adjoint 
operator T* , we have 

(¥>,!) = A 2 , 
(1,^) = 2 • A 3 , 

(<p,tf) = 2-A 2 -exp(-l/(2-A 2 )). 



In particular, 



{(p,l) ■ (1,V) , :i . ,, , 2 , 



A • exp(l/(2 • A )). 



Proof. In the calculations that follows we will use the relations erf(l/(\/2- A)) = \/2/tt, Jq q(x)dx = A 
and ^p*(y)dy = 2-X 2 . 

First the inner product between the eigenfunction and the constant function V. 

((p, 1) = J q(x)dxdy — — ■ j f q(t)dtdxdy 

J[0,1] 2 A Jo< x <y<lJ x 

= [ q(x)dx -\ - ( t ■ (1 - t) ■ q{t)dt 
Jo A J 

= A-A + ^-A 2 -erf(l/(V2-A)) 
v2 

= A 2 . 

Second, the inner product between the adjoint eigenfunction and the constant function 1. We have 
(l,Y>) = / P*{y)dxdy — — - j j p*(u)dudxdy 

J[0,1] 2 A Jo<y<x<\ Jx 

f 1 1 f 1 u 2 

= J p*(y)dy-j-J —-p*{u)du 



2^A ' ' ;! [jl)ihi 



-ro 



The first integral is given by ii = (A - 2 • A 2 - 2 • A 3 ) • exp(l/(2 • A 2 )) - 1/3 + 2 • (A + A 2 + A 3 ). The 
second integral we solve by integration by parts letting /' = ^1 — • \/2 • tt • y • exp anc ^ 

9 = erf (tIa) • Then So ^ gdy = ^0 ~ So f 9 ' dy is given by: 

r.Af „ „ ^ I 1 

h = 



j= ■ A • (2 • A + 2 • A 2 - y 2 ) ■ exp(y 2 /(2 • A 2 )) • erf (y/(V2 • A)) 



J 



" /< 2 

Jo 



A + 2 • A 2 - y 2 )dy. 
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Combining all the terms in the sum I\ + I2 we obtain 2 • A 3 . The third inner product is given by 



(<P,W) = / q{x) ■ p* (y)dxdy 
J[o,i} 2 

— — • f f q(t)dt ■ p* (y)dxdy — — ■ f q(x) ■ f p*(u)dudxdy 

A Jo<x<v<i Jx A Jo<v<x<l Jx 



IO<x<y<Ux a Jo<y<x<i 

1 \ / ,1 x 2 



q(x)dx • / P*(y)dy - T • / t ■ q(t) ■ p*(y)dtdy 

J \Jo J A J0<t<y<l 

2 • A 3 - 2 • A • f (1 - exp(-y 2 /(2 • A 2 ))) • p*(y)dy 
Jo 

-2 • A 3 + 2 • A • f exp(-y 2 /(2 • A 2 )) • p*(y)dy 
Jo 

-2 • A 3 + 4 • A • J (-y + exp(-y 2 /(2 • A 2 )) + yft/ V2 ■ y ■ erf (y/( V2 ■ A))) dy 



2-A 2 -exp(-l/(2-A 2 )). 



□ 



6.2 Asymptotics 



We know that the largest root Ao of the eigenvalue equation (|6.3[) is real, positive and simple, since 
the associated operator T is positivity improving. However, to say a bit more about the eigenvalues 
consider the related equation eii(z) = \/2/tt. 



Since the error function is an increasing function on the real axis, the equation erf (z) = \/2/tt 
has a unique real root. The error function is an odd function hence we know by the strong version 
of the little Picard theorem that the equation erf(z) = y2/vr has infinitely many roots. Moreover, 
the complex roots appear in conjugate pairs. To summarize this discussion we have: The eigenvalue 
equation has a unique real root which is positive and is the largest root. The remaining infinitely 
many roots are all complex and appear in conjugate pairs. 

Numerically, we can approximate the roots of the eigenvalue equation (|6.3p . The unique real root 
is Ao = 0.7839769312 .... The next four largest roots are: 

Ai, 2 = 0.2141426360 . . . ± 0.2085807022 ... -i 
A34 = -0.1677323922 . . . ± 0.2418627350 ... -i 

Furthermore, Proposition 16.21 states that all eigenvalues are simple. Hence from Proposition 16.61 we 
conclude: 

Theorem 6.7. The number of 213 -avoiding permutations satisfies 

^ — (^l)' AS+1 + 0(|Ain 

where Ao = 0.7839769312. . . is the unique real root of the equation erf (l/(\/2 • A)) = \f2/ir, and Ai 
is the next largest root and its modulus is given by \\\\ = 0.298936411 .... 

By considering the next two conjugate roots Ai and A2 we obtain an approximation for the next 
real term in the expansion of o n (213)/n!: 

2 • 1.158597034 . . . • (0.2989364111 . . .) n+1 • cos (-5.593221320 . . . + (n + 1) • 0.7722415374 . . .) . 
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7 123, 231, 312- Avoiding permutations 

Let us consider 123, 231, 312-avoiding permutations. Now the set S is not a singleton, but has cardi- 
nality three, that is, S is given by the set {123,231,312}. The function % is given by 



X(xi,x 2 ,x 3 ) 
That is, the operator T is described by 



1 if X3 < X2 < Xi, 

1 if X2 < X\ < X3, 

1 if X\ < X3 < X2, 

otherwise. 



(Tf)(x,y) 



%f(t,x)dt X0<x<y<l, 
jy f(t, x)dt + £ f(t, x)dt if < y < x < 1. 



7.1 Eigenvalues and eigenfunctions 

For a real number t let {t} denote the fractional part of the real number t, that is, {t} = t — [t\. 
Observe that the fractional part belongs to the interval [0, 1). 

Lemma 7.1. The subspace W of L 2 ([0, l] 2 ) consisting of functions of the form 

f(x,y)=g({x-y}) for (x, y) G [0, l] 2 (7.1) 
is invariant under the operator T. 

Proof. There are two cases to consider. First, if x < y, we have 

T(g({x - y})) = [ V g({t - x})dt = T g(s)ds. 

Jx JO 

Second, if y < x we have 

T(g({x - y})) = [ V g({t - x})dt + f g({t - x})dt 

JO Jx 
rl+y—x rl — x rl+y—x 

= / g(s)ds+ / g(s)ds = / g(s)ds. 

Jl-x Jo Jo 



These two cases can be summarized as the integral from to {y — x}. However the upper bound of 
the integral can be written as 1 — {x — y}, that is, the operator can be expressed as 

l-{x-y} 



ri-^x-y} 

T(g({x - y})) = / g(s)ds. 
Jo 



□ 



Since constant function 1 belongs to W, it is enough to consider eigenvalues and eigenfunctions in 
the subspace W. 
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Proposition 7.2. The eigenvalues of the operator T restricted to the invariant subspace W are given 
by Xj = (— 1)0' —1 )/ 2 ■ 2/(tt ■ j) where j is an odd positive integer and the corresponding eigenf unctions 
are described by 

cp(x,y) = cos({x - y}/\). 



Proof. The defining relations for the eigenfunctions are 



-l-i 

A • g(t) = I g(s)ds. (7.2) 
o 



Observe that A = is not an eigenvalue. Differentiating equation (j7.2|) twice, we obtain 

X-g'(t) = -g(l-t), (7.3) 
A •</'(*) = g'(l-t) = -l/X-g(t). (7.4) 



This second order equation has the general solution A ■ cos(i/A) + B ■ sin(t/A). Setting t = 1 in (|7.2p 
implies that g(l) = and then setting t = in (|7.3|) implies that g'(0) = 0. Hence B = and the 
eigenfunction is cos(t/A). For this function, equation (|T.2[) implies that cos(t/A) = sin((l — t)/X). 
Setting t = we obtain that the eigenvalues satisfy sin(l/A) = 1. They can be parametrized as in the 
proposition. □ 

Note that the eigenvalues tend to as 

Ai > -A 3 > A 5 > -A 7 > • • • . 

Proposition 7.3. The eigenvalues X of the operator T restricted to the subspace W are simple. That 
is, there is no function h(t) such that X ■ h = T(h) + g where g(t) = cos(i/A). 

Proof. The defining relations for the eigenfunctions are 

,i-t 

X-h{t) = / h(s)ds + g(t). (7.5) 



o 



Differentiating this equation twice gives 



X-h'(t) = -h(l-t) + g'(t), (7.6) 
A • h"(t) = h'{l-t)+g"{t). (7.7) 

Substituting (17.61) into (17. 7h and using (17. 4p . we obtain the second order differential equation: 

X 2 -h"{t) + h{t) = 2- X-g"(t). 

The general solution is 

hit) = -4r • t ■ sinit/X) + Ci • cos(t/A) + C 2 • sin(t/A). 
X z 

In the following we use that sin(l/A) = 1 and cos(l/A) = 0. Observe that setting t = 1 in ()7.5[) gives 
X-h(l) = g(l), which implies C 2 = 1/A 2 . Similarly, setting t = in ([75]) gives X-h'(0) = -h(l)+g'(0) 
which implies C2 = 1/(2 • A 2 ), a different value for C 2 . Hence there is no such a function h satisfying 
the integral equation (I7.5p . □ 
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Observe that the function x satisfies the symmetry in Lemma [4 .71 Hence the adjoint eigenfunction 
tp is given by J(p. However, since the value of ip(x, y) only depends on the difference x — y, we have 
that Jip = ip. 

Proposition 7.4. For a non-zero eigenvalue A with eigenfunction ip and eigenfunction ip = if of the 
adjoint operator T* , we have 

(</>,l) = (1,?)=A, 
(<p,ty = 1/2. 



In particular, 



(^i)-M) _ 2 A2 



7.2 Asymptotics 

The above computations show that all eigenvalues of the operator T are simple and given explicitly. 
We thus obtain the following asymptotic expansion for (T n ~ 2 (l),l) = a n (123, 231, 312)/n! as an 
immediate consequence of Theorem 11.11 Propositions 17.21 and 17.41 



" " ' ' 23 ' 231 ' 312) - = 2 • ]T Xj + O (|A 2fc+3 D • 



j odd 



However, when we let k tend to infinity, observe that the right hand side converges. Hence we obtain 
the exact expression 

q n (123,231,312) ^ u N-i)/2-n ( 2 Y 

— 5 — h 

j odd 

for n > 2. However this is the expression for E n _i/(n — 1)! where E n denotes the nth Euler number; 
see equation (|1 .8j) . Thus we conclude that 

Theorem 7.5. The number of 123, 231, 312 -avoiding permutations in & n is given by n ■ E n -% for 
n > 2. 

This result is due to Kitaev and Mansour [14J. Note that the special form of the eigenfunctions and 
the invariant subspace are reflected in their proof that this class of permutations is invariant under 
the shift 7Ti7T2 • ■ • 7T n i — > (tti + 1)(7T2 + !)■•■ (7r n + 1), where the addition is modulo n. 



8 Concluding remarks 

It is straightforward to design a Viennot "pyramid" to compute the number a n of S'-avoiding permu- 
tations. For the original Viennot triangle, see [2D E2]. Let the entry aft"" ,%m of the pyramid be the 
number of permutations in the symmetric group on n elements, avoiding the set S and ending with 
the m entries i±, . . . , i m . Then the entry an''"' 1 " 1 is a sum of entries of the form a > 1 ^y i "' %m ~ 1 . This sum 
is a discrete analogue of the operator T. How far does this analogue between the discrete model and 
the continuous one go? Does the function f n = j 1 ™-"^!) approximate the n-th level of the pyramid? 
More exactly, how well does the integer Un'"' %m compare with n\ ■ f n (ii/n, . . . ,i m /n)7 
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In the case of descent pattern avoidance, can one prove that T restricted to the invariant subspace V 
is compact? We have done so in the case of 123-avoiding permutations. 

Consider the graph G$ of overlapping permutations on the vertex set & m . What is the smallest 
number of edges one has to remove in order to make the graph not strongly connected? Clearly one 
can remove m edges disconnecting the vertex 12 • • • m. Is m the right answer? This would suggest 
that one can remove m — 1 edges without making the directed graph disconnected. 

A more general enumeration problem is as follows. For a function w on 6 m _|_i define the weight of 
a permutation ir = {tt\ti2 ■ ■ ■ vr n ) by the product 

n—m 

Wt(7r) = jj w(II(7rfc,... ,7Tfc +m )). 
k=l 

Now what can be said about the values and asymptotics of the sum 

otiM = ^2 wt ( 7r ) 

TT&6„ 

as n tends to infinity. For instance, if w is a positive function we know by the result of Krem and 
Rutman that 

a„(w) ~ c • A n • nl. 

An operator is called a Volterra operator if its spectral radius is 0. What can be said about Volterra 
operators of the form (II. ID ? More specifically, what can be said about the asymptotic of (T n (l), 1)? 
See Examples II. 13[ fl. 141 and 13.91 for different behaviors. 

In all the cases we have calculated, all the (non-zero) eigenvalues have been simple, see Proposi- 
tions [5T21 [672] and [731 Hence it is natural for us to conjecture: 

Conjecture 8.1. When the graph H$ is ergodic then the associated operator T$ only has simple 
eigenvalues, that is, all the eigenvalues have index 1. 

In the case when associated operator satisfies a symmetry condition as in Lemma 14.71 or in 
Lemma [4. 8 [ does this condition help in proving the conjecture. If Conjecture 18. H is false the asymptotic 
expansion can be obtained by using Theorem 22 in [6] Section VII. 3]. 

In this paper our object is to understand consecutive patterns avoidance. Generalized pattern 
avoidance (also known as vincular pattern avoidance, see Chapter 7 in |13j ) was introduced by Babson 
and Steingrimsson [2]. Is there an analytic approach to obtain asymptotics for these classes of permu- 
tations? Lastly, it would be daring to ask for an analytic proof of the former Stanley- Wilf conjecture, 
recently proved in [17]. 

Acknowledgements 

We are grateful to Judith McDonald for a helpful discussion regarding the Perron-Frobenius Theorem 
and to Hua-Huai (Felix) Chern who helped us to improve Proposition 16.21 The first author thanks 
MIT and Institute for Advanced Study where this paper was completed. The first author was partially 
supported by National Science Foundation grants DMS-0200624 DMS-0902063, DMS-0835373 and 
CCF-0832797, and National Security Agency grant H98230-06- 1-0072. 



34 



References 

[1] S. Avgustinovich and S. Kitaev, On uniquely fe-determined permutations, Discrete Math. 
308 (2008), 1500-1507. 

[2] E. Babson and E. Steingrimsson, Generalized permutation patterns and a classification of 
the Mahonian statistics, Sera. Lothar. Corabin. 44 (2000), Art. B44b, 18 pp. 

[3] M. Bona, Combinatorics of Permutations, Chapman and Hall/CRC Press, 2004. 

[4] A. Claesson, Generalised pattern avoidance, European J. of Combin. 22 (2001), 961-971. 

[5] V. Dotsenko and A. Khoroshkin, Anick-type resolutions and consecutive pattern avoidance, 
arXiv:1002.2761vl [math.CO]. 

[6] N. Dunford and J. T. Schwarz, Linear Operators, Part I: General Theory, New York, Inter- 
science Publishers, 1967. 

[7] R. Ehrenborg and Y. Farjoun, Asymptotics of the Euler number of bipartite graphs, Adv. 
in Appl. Math. 44 (2010), 155-167. 

[8] R. Ehrenborg, M. Levin and M. Readdy, A probabilistic approach to the descent statistic, 
J. Combin. Theory Ser. A. 98 (2002), 150-162. 

[9] S. Elizalde, Asymptotic enumeration of permutations avoiding generalized patterns, Adv. in 
Appl. Math. 36 (2006), 138-155. 

[10] S. Elizalde and M. Noy, Consecutive patterns in permutations, Adv. in Appl. Math. 30 (2003), 
110-125. 

[11] P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009. 
[12] S. Kitaev, Multi-avoidance of generalised patterns, Discrete Math. 260 (2003), 89-100. 
[13] S. Kitaev, Patterns in permutations and words, Springer- Verlag, 2011. 

[14] S. Kitaev and T. Mansour, Simultaneous avoidance of generalized patterns, Ars Combin. 75 
(2005), 267-288. 

[15] D. E. Knuth, The art of computer programming. Vol. 1 Fundamental algorithms., Addison- 
Wesley Publishing Co. Reading, Mass.-London-Don Mills, Ont, 1969. 

[16] M. G. Krein and M. A. Rutman, Linear operators leaving invariant a cone in a Banach space. 
Uspehi Matem. Nauk (N.S.) 3 (1948), no. 1 (23), 3-95. English translation in Amer. Math. Soc. 
Translation 1950 (1950), no. 26, 1-128. 

[17] A. Marcus and G. Tardos, Excluded permutation matrices and the Stanley-Wilf conjecture, 
J. Combin. Theory Ser. A. 107 (2004), 153-160. 

[18] A. Mendes and J. Remmel, Permutations and words counted by consecutive patterns, Adv. in 
Appl. Math. 37 (2006), 443-480. 

[19] R. Simion and F. Schmidt, Restricted permutations, European J. Combin. 6 (1985), 383-406. 



35 



[20] R. P. Stanley, Enumerative Combinatorics, Vol. 1, Wadsworth & Brooks/Cole, 1986. 

[21] G. Viennot, Permutations ayant une forme donnee, Discrete Math. 26 (1979), 279-284. 

[22] G. Viennot, Equidistribution des permutations ayant une forme donnee selon les avances et 
coavances, J. Combin. Theory Ser. A. 31 (1981), 43-55. 

[23] R. Warlimont, Permutations avoiding consecutive patterns, Ann. Univ. Sci. Budapest. Sect. 
Comput. 22 (2003), 373-393. 

R. Ehrenborg, and P. Perry, Department of Mathematics, University of Kentucky, Lexington, KY 
40506-0027, USA, { jrge,perry}@ms . uky.edu 

S. Kitaev, Department of Computer and Information Sciences, University of Strathclyde, Livingstone 
Tower, 26 Richmond Street, Glasgow Gl 1XH, United Kingdom, sergey.kitaev@cis.strath.ac.uk 



36 



